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SUMMARY 


As  a  part  of  a  larger  research  program  to  bring  about  acceptability  of 
fiber-reinforced  plastics  for  prhnary  load-bearing  members  of  Army  aircraft 
structures,  a  series  of  mechanical  tests  was  conducted  on  fiher  glass-epcxy 
laminates.  The  loadings  used  were  uniaxial  tension;  biaxial  tension  with 
principal-stress  ratios  of  1:2,  1:1,  and  2:1;  and  pure  shear.  The  most 
unique  tests  were  those  involving  the  biaxial  loadings.  These  tests  we  e 
carried  out  on  a  specially  designed  test  fixture  using  a  cable  pulley  sys¬ 
tem  in  conjunction  with  a  conventional  universal  testing  machine.  Data 
were  obtained  on  the  limit  strength  and  apparent  ultimate  strength  of  fiat 
laminates  subjected  to  biaxial  loading.  Under  biaxial  loading,  the 
parallel-ply  lamination  arrangement  was  found  to  be  > lightly  superior  to 
the  cross-ply  and  quasi- isotropic  lamination  arrangements.  The  shear  tests 
were  accomplished  on  thin-walled  tubular  specimens  subjected  to  pure  tor¬ 
sion. 
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FOREWORD 


This  report  was  prepared  by  the  University  of  Oklahoma  Research  Institute 
(OURI)  under  Phase  II  of  U.S.  Army  Aviation  Materiel  Laboratories  (USAAVLABS) 
Contract  DAAJQ2-67-C-Q11I  for  research  conducted  during  the  period  from  July 
1,  1967  to  July  31,  1968.  Work  under  this  contract  has  also  included  sandwich- 
shell  vibrational  investigations  (Phase  I),  reported  separately  in  Technical 
Report  68-85,  titled  Vibration  Evaluation  of  Sandwich  Conical  Shells  with 
Fiber-Reinforced  Composite  Facings! 

The  research  effort  is  a  continuation  of  the  research  accomplished  under 
three  previous  contracts  and  reported  in  the  following  USAAVLABS  reports: 

1.  Technical  Report  67-65,  Fabrication  and  Full  Scale  Structural 
Evaluation  of  Sandwich  Shells  of  Revolution  Composed  of  Fiber 
tiTass  Reinforced  Piastre  Facings  ail'd  Honeycomb'  Cores'. 

2.  Technical  Report  65-66,  The  Effect  of  Resin  Content  and  Voids  on 
the  Strength  of  Fiberglass-Reinforced  Plastics'  for  Airframe  Use. 

3.  Technical  Report  65-60,  Dynamic  Elastic  Damping,  and  Fatigue 
Characteristics  of  Fiberglass -Reinforced  Sandwich  Structure. 

4.  Technical  Report  65-15,  Strength  Properties  and  Relationships 
Associated  with  Various  Types  of  Fiberglass -Reinforced-Pacing 
Sandwich  Structured 

5.  Technical  Report  64-37,  Research  in  the  Field  of  Fiberglass- 
Reinforced  Sandwich  Structure  for  Airframe  Use. 

The  present  report  was  written  by  Dr.  Charles  W.  Bert,  project  director  and 
Professor  of  Aerospace  and  Mechanical  Engineering;  Mr.  Byron  L.  Mayberry, 
research  engineer;  and  Dr.  John  D.  Ray,  project  director. 

Special  acknowledgment  is  made  to  Mr.  Frederick  Lehmann  in  connection  with 
fabricating  the  specimens,  and  to  Mr.  Jack  Holland  in  fabricating  the 
special  test  equipment  and  in  performing  most  of  the  experiments. 
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1.  INTRODUCTION 


Under  the  sponsorship  of  the  U.S.  Army  Aviation  Materiel  Laboratories 
(USAAVLABS) ,  the  University  of  Oklahoma  Research  Institute  (QURI)  has  been 
conducting  A  major  research  program  for  the  past  five  years  in  the  area  of 
fiber  glass  reinforced  plastics  (FRP)  suitable  for  primary  structures  of 
U.S.  Army  aircraft.  Tasks  whicn  have  been  completed  include  the  evaluation 
of 

(1)  Several  sandwich  fabrication  techniques  (Reference  1) . 

(2)  Effect  of  primary  fabrication  variables  on  conventional  strength  proper¬ 
ties  of  laminates  (References  1,  2,  and  3). 

(3)  Several  honeycomb -core  sandwich  configurations  as  small  panels  sized  to 
fail  in  face  rupture  and  various  buckling  modes  (References  1,  4,  and  5). 

(4)  Fatigue  properties  of  sandwich  material  as  beams  with  FRP  facings  and 
honeycomb  core  (References  6  and  7) . 

(5)  Dynamic  stiffness  and  damping  of  sandwich  material  in  the  form  of  beams 
(References  6  and  8) . 

(6)  Fabrication  and  static  buckling  of  full-size  sandwich  shell  structures 
in  the  form  of  cylindrically  curved  panels,  complete  cylinders,  and 
truncated  cones  (References  9,  10,  11,  and  12). 

(7)  Vibrational  frequencies,  modes,  and  damping  of  full-size,  sandwich, 
truncated  conical  shells  (References  13  and  14). 

The  objectives  of  the  present  research  were  (1)  to  develop  equipment  and 
procedures  for  detexmining  the  mechanical  behavior  of  laminates,  (2)  to 
evaluate  the  mechanical  behavior  of  fiber  glass  reinforced  plastic  (FRP) 
laminates  under  uniaxial  tension,  biaxial  tension  (1:2,  1:1,  and  2:1  prin¬ 
cipal-load  ratios),  and  pure  shear,  aid  (3)  to  relate  this  behavior  to 
laminate  mechanics  analysis  methods. 

For  simplicity,  most  previous  evaluations  of  the  mechanical  behavior  of 
composite  materials  were  limited  to  simple  loadings,  i.e. ,  uniaxial  tension 
and  compression,  flexure,  and  pure  shear  (References  15  and  16).  In  sane 
of  these  investigations,  the  very  significant  effects  of  angular  orientation 
of  the  loading  with  respect  to  the  major  material- symmetry  axis  of  the  com¬ 
posite  material  have  been  studied  (References  17  and  18).  However,  in 
many  regions  of  practical  aircraft  structures,  the  material  is  subjected  to 
biaxial  loadings,  i.e.,  loads  in  two  directions. 

There  is  a  large  amount  of  biaxial- load  test  data  for  metallic  alloys  used 
in  aircraft  structures.  However,  in  contrast,  only  a  very  limited  amount 
of  such  data  is  available  for  composite  materials. 
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In  the  present  program,  an  original  cable-and-pulley  system  reacted  by  a 
rigid  frame  was  used  to  load  the  specimen  biaxially.  With  this  simple 
arrangement,  described  in  Section  5a,  the  specimen  is  loaded  biaxially, 
although  the  external  load  is  applied  in  one  direction  only. 

The  material  used  was  181-style  E-glass  cloth  impregnated  with  Epon-828 
epoxy  and  activated  by  curing  agent  Z.  In  Table  I  are  tabulated  the 
lamination  arrangements  and  loading  orientations  used. 


TABLE  I.  LAMINATION  ARRANGEMENTS  AND  LOADINGS 

Lamination 

Type  of 

No.  of 

Angular  Orientation 

Arrangement 

Loading 

Plies 

of  Loading,  deg 

Parallel  ply 

Biaxial  1:1 

4 

0  8  45 

Biaxial  1:2 

4 

0  8  90* 

Shear 

4 

0  6  45 

Uniaxial  tension 

4 

0  6  90 

Cross  ply 

Biaxial  1:1 

4 

0 

(balanced) 

Biaxial  1:2 

4 

0** 

Shear 

4 

0  8  45 

Uniaxial  tension 

4 

0** 

Quasi- isotropic  Biaxial  1:1 

6 

0  8  45 

(0  -60  -120") 

Biaxial  1:2 

6 

0  8  90* 

Shear 

6 

0  8  45 

Uniaxial  tension 

6 

0  8  90 

*A  1:2  biaxial 

loading  at  90  deg  is 

sometimes 

designated  at  a  2:1  biaxial 

loading  at  0  deg. 

**In  a  balanced  cross-ply  laminate. 

there  is  no  difference  between  the  effect 

of  a  loading  at  90  deg  and  that  of  the  same 

loading  at  0  deg.  Thus,  tests 

at  only  0  deg  needed  to  be  conducted. 

Five  tests  were  carried  out  for  each  combination  tabulated  above.  Thus,  there 
were  50  biaxial  tests,  30  shear  tests,  and  25  uniaxial  tension  tests,  or  a 
total  of  105  tests. 
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THEORY  OF  THE  IN- PLANE  STIFFNESSES  OF  COMPOSITE  PLATES 
iSYNflET'RIOAIXY  LMMaTEH  OF'  oRTHoTROPIC  UVERS 


The  general  theory  presented  here  starts  in  Section  2a  with  development 
from  first  principles  of  the  general  three-dimensional  theory  of  aniso¬ 
tropic  elasticity.  Then  the  general  theory  is  reduced  to  the  practical 
case  of  a  thin  orthotropic  layer,  such  as  one  layer  in  an  FRP  composite. 

In  Section  2b  the  theory  is  presented  for  a  plate  loaded  in  its  plane  only 
and  consisting  of  an  arbitrary  number  of  sets  of  orthotropic  plies.  Each 
set  consists  of  an  arbitrary  number  of  pairs  of  plies,  and  each  pair  has 
one  ply  located  a  certain  distance  above  the  midplane  of  the  laminated 
plate  and  the  other  ply  located  the  same  distance  below  the  midplane.  All 
plies  in  a  given  set  have  the  same  orientation. 

a.  A  Single  Orthotropic  Layer 

The  currently  popular  idea  of  considering  a  composite  material,  con¬ 
sisting  of  oriented  fibers  imbedded  in  a  matrix,  to  behave  macrosco- 
pically  as  a  homogeneous,  anisotropic  material  is  used  hereT  However, 
this  is  not  new,  since  it  was  used  as  early  as  1914  by  Huber  (Reference 
19)  in  connection  with  reinforced  concrete  slabs. 

The  anisotropic  elastic  coefficients  may  be  determined  either  experi¬ 
mentally  or  analytically.  In  the  past  decade,  literally  dozens  of 
papers  have  been  published  on  theoretical  micro-mechanics  analyses  for 
predicting  the  anisotropic  elastic  coefficients  of  single  layers 
reinforced  with  unidirectional  fibers  (such  as  obtained  in  che  filament¬ 
winding  process) .  Many  of  these  have  been  verified  experimentally  to 
be  sufficiently  accurate  for  engineering  design.  In  contrast,  the 
authors  know  of  no  theoretical  analysis  for  predicting  the  stiffness  of 
a  single  layer  of  a  woven-fabric  reinforced  composite.  However,  it  is 
well  known  that  due  to  the  crimp  of  the  fibers  in  a  woven-cloth  rein¬ 
forced  composite,  such  composites  have  lower  stiffnesses  than  unidirec- 
tionally  reinforced  composites  with  the  same  fiber  content.  Therefore, 
in  the  present  work,  the  anisotropic  elastic  coefficients  are  determined 
purely  by  means  of  experiments  (see  later  sections) . 

Anisotropic  materials  have  material  properties  which  depend  upon  the 
directional  orientation.  The  most  general  stress-strain  relation  for  a 
linear  e^stic  material  is  known  as  generalized  Hooke's  law.  For  a 
general  three-dimensional  anisotropic  body,  generalized  Hooke's  law  can 
be  expresseo  in  a  Cartesian  coordinate  system  (x,  y)  as  follows 
(Reference  20) : 
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The  oTg  are  the  stress  components ,  the  e|g  are  the  strain  components, 
and  the  Cj4  are  the  Cauchy  stiffness  coefficients.  The  subscripts  on 
the  stresses  are  determined  as  follows:  the  first  subscript  denotes 
the  direction  of  the  noitnal  to  the  plane  on  which  the  stress  acts ,  and 
the  second  subscript  denotes  the  direction  in  which  the  stress  component 
acts.  For  example,  a--  denotes  a  stress  acting  on  a  plane  perpendicular 
to  the  x  axis  and  in  tne  direction  of  the  y  axis.  A  similar  notation 
holds  for  the  strain  components;  i.e.,  the  first  subscript  denotes  the 
direction  of  the  normal  to  the  plane  being  strained,  and  the  second  sub¬ 
script  denotes  the  direction  in  which  the  strain  is  taking  place. 

Of  the  six  stress  components,  the  three  having  both  subscripts  the  same 
(xx,  yy,  zz)  are  oriented  in  a  direction  normal  to  the  surface  on  which 
they  are  acting  andthus  are  cal  led  normal  stresses.  The  other  three 
stress  components  (yz,  zx~  xy)  are  oriented  in  a  direction  tangential 
to  the  surface  on  which  they  act  and  are  known  as  shear  stresses. 

It  is  noted  that  the  subscript  notation  used  in  Equation  (1)  for  the 
Cauchy  coefficients  is  dependent  only  on  the  position  of  the  coefficient 
in  the  array.  The  particular  order  in  which  the  elements  appear  in  the 
stress  and  strain  colunn  matrixes  is  fairly  well  standardized,  having 
been  used  many  years  ago  by  A.E.H.  Love;  Reference  20  is  the  last  edition 
of  his  book. 


Equation  (1)  can  be  rewritten  more  compactly  in  matrix  notation  as 
follows: 


COj]  =  [CyHtjJ  (2) 

where  [ajj  and  [e-j]  are  the  column  vectors  representing  the  respective 
stress  and  strain  components,  C^j  is  the  square  matrix  of  Cauchy  stiff¬ 
ness  coefficients,  and  the  subscripts  of  the  stress  and  strain  components 
have  been  redesignated  as  follows: 

xx  -+  1,  yy  +2,  zz  -*■  3, 

yz  ■*  4,  zx  -*•  5,  xy  -*■  6 
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An  alternate,  but  more  complicated,  notation  for  the  stiffness  coeffi¬ 
cients  is  Cartesian  tensor  notation.  However,  this  requires  four  sub¬ 
scripts  on  each  coefficient*,  and  thus  the  simpler  double- subscript 
matrix  notation  is  used  here. 


Assuming  that  there  are  no  electromechanical  effects  and  that  the  defor¬ 
mation  takes  place  isothermally,  it  can  be  shown  by  use  of  thermodynamics 
(Reference  22)  that  Castigliano's  theorems  must  hold  for  infinitesimal 
displacements  of  a  linear  elastic  material.**  These  theorems  can  be 
expressed  mathematically  as  follows : 


5U/3E fg  *  °fg 

3U/80fg  = 


(3) 


where  U  is  the  potential  strain  energy  per  unit  volume, 

u  =  I  /  afg  d£fg 

Equations  (3)  require  that  the  matrix  of  Cauchy  coefficients  be  symme¬ 
tric,  i.e., 


(4) 


A  single  layer  of  any  known  composite  material  has  certain  planes  of 
elastic  symmetry.  In  a  layer  consisting  of  unidirectionally  oriented 
fibers,  the  material- symmetry  planes  are  the  plane  normal  to  the  fibers 
and  any  two  orthogonal  planes  parallel  to  the  fibers.  In  a  perpendicu¬ 
larly  woven- cloth  layer  such  as  used  in  the  experiments  reported  here, 
the  material- symmetry  planes  are  the  two  planes  normal  to  the  two  fiber 
orientations  and  the  plane  of  the  layer  itself.  Such  a  system  of  three 
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*The  general  stiffness-coefficient  tensor  would  have  (3)  =  81  elements. 

However,  neglecting  microstructural  effects  (Reference  21)  as  is  customary 
in  classical  theory  of  elasticity,  equilibrium  requires  that  the  stress 
tensor  be  symmetric  (a  r  =  Or  )  and  infinitesimal-displacement  kinematics 
requires  that  the  strain  tenser  be  symmetric  (e  r  =  eJ .  Thus,  the 
coefficient  tensor  would  have  only  36  independent  elements,  in  agreement 
with  Equation  (1) . 

**In  the  case  of  geometrical  nonlinearity  (large  displacements)  or  material 
nonlinearity  (nonlinear  stress-strain  relation),  the  complimentary  energy 
must  be  used  instead  of  the  potential  energy  (Reference  23) . 
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orthogonal  planes  of  material  symmetry  is  said  to  be  rectangularly 
orthotropic.*  The  intersections  of  the  three  orthogonal  planes  of 
symmetry  make  three  orthogonal  axes  referred  to  as  the  material  -  symmetry 
directions .  If  the  Cartesian  coordinate  axes  are  selected  to  coincide 
with  the  material -symmetry  directions,  24  of  the  36  coefficients  appear¬ 
ing  in  the  stiffness -coefficient  array  in  Equation  (1)  disappear.  Incor¬ 
porating  this  fact,  as  well  as  the  symmetry  relation,  Equation  (4),  the 
three-dimensional  array  becomes: 
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To  reduce  relation  (5)  to  a  two-dimensional  case,  a  state  of  generalized 
plane  stress  is  assuned;  i.e.,  the  layer  is  assumed  to  be  sufficiently 
thin  compared  to  its  other  dimensions  that  all  three  components  in  the 
thickness  direction  are  negligible.  Thus, 


o  =  0  =  CT1  e  —  +  C,~  £ —  +  C,7  e — 
31  xx  32  yy  33  zz 
zz  ' ' 
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(6) 


°zx  ^  ^55  ezx 

The  last  two  of  Equations  (6)  imply  that  the  thickness-direction  shear 
strains  e--  and  £--  are  zero. 

Although  the  thickness -direction  normal  stress  a^;  is  zero,  the  thickness- 
direction  noimal  strain  e-~  is  not  zero  in  general.  However,  e~  may  be 
expressed  in  terms  of  e^x  and  £yy  by  means  of  the  first  of  Equations  (6). 
For  a  thin  orthotropic  layer  with  the  coordinate  axes  x,  y  coinciding 
with  the  material -symmetry  directions,  the  stress-strain  relations  are 


*  More  complicated  classes  of  elastic  symmetry  than  orthotropic  are  found  in 
single  crystals  of  crystalline  materials  (Reference  24),  and  tree  trunks 
are  examples  of  cylindrically  orthotropic  material  (Reference  25).  Also, 
in  micromechanics  analyses  of  fiber- reinforced  composites,  hexagonal  sym¬ 
metry  is  sometimes  assuned  (Reference  26).  However,  the  macroscopic 
behavior  of  a  thin  layer  of  fiber -reinforced  composite  material  is  repre¬ 
sented  adequately  as  a  rectangularly  orthotropic  material. 
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as  follows: 
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where  the  reduced  stiffness  coefficients  are  related  to  the  three-dimen¬ 
sional  Cauchy  coefficients  as  follows: 

^11  =  C11  ‘  ^C13^C33^ 

^12  =  C12  ‘  ^C13C23^C33^  (8) 

^66  =  C66 


It  is  noted  that  in  the  array  of  five  non-zero  reduced-stiffness 
coefficients  appearing  ^n  Equation  (7),  only  four  are  independent.* 


Rather  than  using  the  stiffness  coefficients,  most  design  and  structural 
engineers  use  the  so-called  engineering  moduli:  Young's  elastic  moduli 
E-q  and  E22,  Poisson’s  ratios  and  anc*  shear  modulus  G.  They 
are  most  easily  defined  by  writing  the  strains  in  terms  of  the  stresses 
as  follows: 
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In  view  of  the  symmetry  resulting  from  the  general  reciprocal  relation. 
Equation  (4) , 


v21^E22  v12^E11 

Then  one  of  the  Poisson's  ratios  can  be  expressed  in  terms  of  the  other 
and  the  two  Young's  moduli.  Thus,  the  four  independent  engineering 
elastic  moduli  for  a  thin  orthotropic  layer  are  E^,  E22,  v^,  and  G. 

*It  is  often  assumed  that  fiber-reinforced  composites  (especially  filament- 
wound  ones)  behave  isotropically  in  the  plane  normal  to  the  fibers.  Thus, 
they  are  said  to  be  "transversal ly  isotropic",  although  a  better  tern  would 
be  "mono tropic"  (Reference  27) .  However,  since  no  thickness -direction 
quantities  enter  into  the  stress-strain  relations  for  a  thin  member,  there 
is  no  difference  between  the  elastic  coefficients  of  such  a  material  and 
an  orthotropic  one. 
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Ordinary  isotropic  materials  are  a  special  case  of  orthotropic  materials. 
However,  it  is  emphasized  that  there  are  two  necessary  and  sufficient 
conditions  for  a  material  to  be  isotropic: 

E22  =  E11  =  E 

(This  implies  that  =  v) 

E/G  =  2(1  +  v)  (12) 

Unidirectional ly  reinforced  composites  usually  have  E~~  <<  E.,.  In 
contrast,  a  bidirectionally  reinforced  layer,  such  as  a  single  layer  of 
woven-cloth-reinforced  composite,  usually  has  E_?  =  E-..  However,  this 
does  not  necessarily  mean  that  Equation  (12)  is^satisfied.  In  fact,  for 
most  cloth-reinforced  composites,  E/G  >  2(1  +  v) . 


In  experimental  stress  analyses,  in  which- strains  rather  than  stresses 
are  measured  experimentally,  it  is  necessary  to  invert  Equations  (9) , 
with  the  following  result: 
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where 
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1  •  v12v21 


(14) 


Direct  comparison  of  the  coefficients  appearing  in  Equations  (7)  and 
(14)  yields  the  following  relations  among  the  reduced-stiffness  coeffi¬ 
cients  and  the  engineering  moduli: 


^11  =  E11^X  *  ^12  =  v21E11'/X  S!  v12E22//X  ’  ^22  “  E22^X 

%t-G 


(15) 


Sometimes  it  is  necessary,  even  in  a  single-orthotropic  layer  structure, 
to  refer  the  elastic  coefficients  to  axes  which  do  not  coincide  with  the 
material-symmetry  axes.  An  example  of  this  is  a  rectangular  plate  con¬ 
sisting  of  a  single  orthotropic  layer  having  its  major  material-symmetry 
axis  oriented  at  an  acute  angle  with  the  longest  side  of  the  plate.  Then  the 
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stress  tensor  and  the  strain  tensor*  transform  as  a  second-rank  tensor, 
and  the  Mohr  stress  and  strain  circles  may  he  used  in  the  two-dimensional 
case.  The  transformed  stress-strain  relations,  related  to  new  orthogonal 
axes  x,  y  which  do  not  coincide  with,  the  material -symmetry  axes,  are 
written  as  follows: 
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where  m  e  cos  0^,  n  s  sin  0^.  It  is  noted  that  since  the  layer  is 
orthotropic,  and  do  not  exist.  However,  in  general,  this  does 
not  imply  tharQj^  and^Q^  are  zero  also. 


One  of  the  most  distinct  differences  between  the  elastic  behavior  of  an 
orthotropic  material  and  that  of  one  which  is  isotropic  is  in  connection 
with  the  directions  of  principal  stresses  and  of  principal  strains.  In 
isotropic  materials,  it  is  well  known  that  the  directions  of  principal 
stresses  and  of  principal  strains  coincide.  However,  as  pointed  out 
very  clearly  by  Greszczuk  (Reference  28),  this  is  generally  not  the  case 
in  orthotropic  materials.  In  practical  numerical  examples  he  cited,  the 
two  principal  directions  differed  by  as  much  as  15  to  30  degrees. 


*The  ordinary  stress  canponents  foim  the  stress  tensor.  However,  in  order 
for  the  strain  components  to  make  up  a  tensor,  the  ordinary  shear  strains 
(such  as  used  here)  must  be  multiplied  by  a  factor  of  one  half.  This 
accounts  for  the  necessity  of  introducing  a  factor  of  one  half  to  shear 
strains  on  the  Mohr’s  strain  circle. 
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b.  General  Symmetrically  Laminated  Plates 

The  theory  of  stiffness  of  laminated  plates  presented  in  this  section 
generally  follows  the  work  of  Tsai  (References  17  and  29),  except  that 
in  Section  2b,  provision  is  made  for  symmetrically  laminated  sets  of 
plies.  All  plies  of  a  given  set  are  assumed  to  have  the  same  properties 
and  orientation.  By  a  symmetrically  laminated  set,  we  mean  one  which 
has  corresponding  pairs  of  plies  located  above  and  below  the  midplane 
of  the  plate.  Then  there  can  be  no  coupling  between  bending/ twisting 
effects  and  in-plane  (stretching/plane  shear)  effects.*  Thus,  when  a 
symmetrically  laminated  plate  is  subject  to  in-plane  loading,  no  bend¬ 
ing  effects  are  induced. 

Generally  a  laminated  plate  has  different  stresses  in  different  sets  of 
plies.  Thus,  one  cannot  speak  of  a  single  stress  at  a  given  plate  cross 
section;  it  is  necessary  to  speak  of  N,  the  load  resultant  expressed  in 
force  per  "length  of  run"  (width) ,  which  is  the  combined  effect  of  the 
stresses  in  the  various  sets  of  plies.  In  general, 


h/2 
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(18) 


where  h  is  the  total  plate  thickness,  and  z  is  the  distance  from  the 
midplane  of  the  plate. 

The  in-plane  stiffnesses  of  a  laminated  plate  are  defined  as  follows: 
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where,  strictly  speaking,  the  e  refers  to  the  strain  components  in  the 
midplane.  However,  since  the  only  loading  considered  here  is  in-plane 
loading  and  since  there  are  no  bending  or  twisting  effects  induced  by 
coupling,  within  the  framework  of  thin-plate  theory,  there  can  be  no 
strain  gradients  through  the  thickness.  In  other  words,  within  the 
assumptions  of  the  present  theory,  the  in-plane  strain  components  can¬ 
not  vary  through  the  thickness. 


In  view  of  the  Equations  (18)  and  (19)  and  the  subsequent  discussion,  it 
is  clear  Jiat  the  general  equation  for  calculating  the  composite  stiff¬ 
ness  coefficients  is  as  follows: 


*)For  treatment  of  such  coupling  effects,  the  reader  is  referred  to  the  works 
of  Stavsky  (Reference  30)  and  Tsai  (Reference  17). 
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(20) 
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where  n  refers  to  a  typical  individual  ply,  N  is  the  total  number  of 
plies,  and  Qjj(n)  is  the  reduced  elastic  coefficient  of  the  nth  ply. 

Assuming  that  the  elastic  coefficients  of  each  individual  ply  are  uni¬ 
form  through  the  thickness  of  the  ply  and  that  there  are  K  sets  of  plies 
(k  being  a  typical  set),  with  all  plies  in  a  given  set  having  the  same 
properties  oriented  in  the  same  direction,  Equation  (20)  becomes 


K  N/K 
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(21) 


If  it  is  further  assumed  that  all  plies  have  the  same  thickness  and 
that  all  K  sets  have  the  same  number  of  plies,  Equation  (21)  can  be 
simplified  still  further  as  follows: 


Ajj  =  (h/K) 
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(22) 


Sometimes  it  is  convenient  to  express  the  composite  coefficients  in  the 
form  of  engineering  properties  defined  as  follows : 

®11  =  t1  *  v12v21^All/h) 

^22  =  ^  '  v12v21^A22^h^ 

*  -  V*  (23) 

v12  =  A12/A22 


In  all  of  the  theory  presented  here,  the  interply  shear  strain  energy 
has  been  neglected.  This  strain  energy  is  due  to  relative  rotation 
between  two  adjacent  plies,  and  it  has  been  treated  by  Clark  and  Chamis 
(References  16  and  31) .  However,  for  the  lamination  arrangements  used 
in  the  present  research,  this  effect  would  be  relatively  small. 

Tsai  (Reference  18)  has  shown  the  importance  of  the  residual  stresses, 
caused  by  thermoelastic  effects  produced  by  the  difference  between  the 
laminating  temperature  and  the  operating  temperature,  in  determining 
failure  of  laminates.  However,  since  all  of  the  laminates  used  in  the 
present  program  were  made  under  closely  controlled  conditions  and 
tested  at  temperatures  which  varied  very  little,  this  effect  would  be 
the  .same  (i.e.,  a  constant  factor)  in  all  of  the  tests  covered.  Thus, 
in  the  theory  presented  here,  thermoelastic  effects  are  emitted  for 
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brevity.  More  detailed  information  on  this  topic  is  given  in  References 
16  and  32. 

c.  Parallel-Ply  Laminates 

Any  parallel-ply  laminate  consisting  of  identical  layers  is  obviously 
symmetrically  laminated.  Then  K  =  1  and  Equation  (22)  becomes 


(24) 


It  is  convenient  to  select  the  composite  reference  axes  to  coincide 
with  the  material-synmetry  axes.  Then  0  =  0  and  Equation  (24)  can  be 
rewritten  as  follows: 


Specifically, 


A12  =  h  Q12  =  h  v12  E22/X 
A22  =  h  ^22  =  h  E22^X 
A66  =  h  Q66  =  h  G;  A16  =  A26  =  0 

Obviously, 

E11  =  Ell»  E22  =  E22»  E  =  G»  v12  =  v12 


(25) 


(26) 


(27) 


It  is  noted  that  a  parallel-ply  laminate  consisting  of  identical  ortho¬ 
tropic  layers  behaves  orthotropically  also. 


d.  Cross-Ply  Laminates 

The  smallest  number  of  plies  with  which  it  is  possible  to  have  a  symme¬ 
trically  laminated  cross-ply  laminate  is  four,  as  shown  schematically 
in  Figure  1. 


The  cross-ply  laminate  consists  of  two  sets  of  layers,  i.e.,  K  =  2. 

The  layers  aligned  in  the  same  direction  as  the  outer  layers  (in  the 
case  of  a  four-ply  laminate,  these  are  the  outer  layers  only)  are  de¬ 
signated  as  set  1,  and  the  other  layers  (i.e.,  the  perpendicular  layers) 
are  designated  as  set  2.  Then  Equation  (22)  becomes 

Ay  -  (h/2)  [QI,™  ♦  Q},®]  (28) 
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Figure  1.  Schematic  Diagram  of  Symmetrically  Laminated  Cross-Ply  Composite. 


For  convenience,  the  direction  of  major  material  symmetry  of  set  1  is 
designated  as  direction  1  of  the  composite.  Then 


Specifically, 


Q!  =  Q 

Q!  .(2)  =  Q.  . 
yij  yr) 


(29) 


Then 


(h/2)(Qn  +  Q22)  =  (h/2)(En  +  E22)/(l  -  v12  v21) 

(h/2) (Q12  +  Q21)  =  h  Q12  =  h  v12E22/(1  -  v12  v21) 

(h/2)  (Q22  +  Qn)  =  (h/2)(E22  +  Eu)/(1  -  v12  v21)  =  AR  (30) 

(h/2)  (Q66  +  Q^)  =  h  =  hG 


*11  =  ^/2^EH  +  **22*  =  E22»  G  G 
v12  =  v12  (2E22)/(E11  +  E22* 


It  is  noted  that  a  symmetrically  laminated  cross-ply  laminate  has  iden¬ 
tical  Young's  moduli  (ft, ,  =  ,  but  it  nevertheless  behaves  orthotro- 
pically  (rather  than  isotropically)  because  ft  f  E*/  [2(1  +  v)j. 

Care  must  be  exercised  in  interpreting  strain-gage  data  obtained  on 
cross-ply  laminates.  This  is  shown  by  the  following  analysis.  For  set  1, 
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A  similar  group  of  relationships  holds  for  set  2. 


(32) 


To  determine  the  actual  biaxial- load  ratio  attained  in  the  test  section, 
another  set  of  relations  is  necessary.  For  gages  aligned  with  the 
material-synmetry  axes, 


N11  A11  ell  +  A12  e22 

N22  =  A12  £11  +  A22  e22 

Thus 

N22^N11  =  ^A12  ell  +  A22  g22^A11  g11  +  A12  e22^ 
or 

N22/^N11  =  ^A12  +  A22  ^22^11^/^11  +  A12  ^22^^11^ ^ 
or 

N22^N11  =  ^v12  +  ^e22^ell^  1 1  +  V12  ^e22^Gll^  ^  (33) 


e.  Quasi- Isotropic  Laminates 

Werren  and  Norris  (Reference  33)  have  shown  mathematically  that  it  is 
possible  to  orient  the  individual  layers  (or  sets  of  layers)  in  multiple- 
layer  laminates  in  such  a  way  that  the  resulting  in-plane  elastic 
behavior  is  isotropic,  i.e.,  has  elastic  coefficients  which  are  indepen¬ 
dent  of  orientation  in  the  plane.  The  conditions  which  must  be  met  by 
the  laminating  arrangement  to  achieve  this  isotropic  behavior  are  as 
follows:  (1)  there  must  be  three  or  more  sets  of  layers,  (2)  the  indi¬ 
vidual  layers  must  have  identical  thicknesses  and  identical  orthotropic 
elastic  coefficients  with  respect  to  their  material -symmetry  axes,  and 
(3)  a  typical  set  of  layers  k  (ranging  from  one  to  K)  must  be  oriented 
at  the  following  angle  with  respect  to  a  reference  direction: 

ek  =  tt  (k  -  1)/K  (34) 

In  the  light  of  the  coupled  theory  (References  17  and  30) ,  the  Werren- 
Norris  design  would  be  isotropic  with  respect  to  in-plane  stiffness, 


14 


but  not  with  respect  to  bending  stiffness  and  bending/in-plane  coupling. 
Thus,  the  Werren-Norris  design  is  called  quasi-isotropic  here. 

To  apply  Equation  (22)  with  the  Q!  given  by  Equation  (17)  and  the 
0i  given  hy  Equation  (34),  certain  sums  presented  by  Werren  and  Norris 
(Reference  33)  are  needed: 

K  K 

y  cos4  [n(k  -  1)/K]  =  l  sin4  [ir(k  -  1)/K]  =  3K/8 
k=l  k=l 


K 

l  sin2  [tt  (k  -  1)/K]  cos2  [tt (k  -  1)/K]  =  K/8 
k=l 

K 

l  cos3  [tt (k  -  1)/K]  cos  [tt (k  -  1)/K] 
k=l 

K 

=  l  cos  [tt (k  -  1)/K]  sin3  [tt (k  -  1)/K]  =  0 
k=l 


(35) 


The  elements  of  the  in-plane  stiffness  matrix  are  given  by  the  following 
expressions: 

An/h  =  A22/h  =  (1/8)[3(Q11  +  Q22)  +  2  (Q12  +  ^66^ 

A12/h  =  (1/8)  (6  Q12  ♦  Qn  ♦  Q22  -  4  Q66) 

A66/h-  (1/8J(4  -  QX1  -  Q22  -  2  Q12J  (36) 

Al6/h  =  A26/h  =  0 

It  is  noted  that  since  A15  =  A26  =  0,  the  in-plane  stiffness  exhibits 
at  least  orthotropic  symmetry.  Furthermore,  since  A^  =  Aj2  and  A^  = 
(1/2)  (A^i  -  A^),  the  in-plane  stiffness  has  complete  symmetry,  i.e., 
it  is  isotropic.  It  is  important  to  notice  that  the  Aji  are  unaffected 
by  the  order  of  lamination;  i.e.,  the  (n  +  1)  layer  does  not  have  to  be 
next  to  the  n  layer,  etc.  Also,  the  Ajj  are  independent  of  the  number 
of  layers  N. 

Converting  the  above  expressions,  Equations  (36),  to  engineering  moduli 
yields  the  following  results: 

v  =  (E-j^  +  E22  +  6v2jE^  -  4XG)/(3En  +  3E22  +  Zv^E^  +  4iG) 

K  =  (1/8)  (1  -  v'2)x"1(3En  +  3E22  +  2v21En  +  4XG)  (37) 

tr  =  f/[2(i  +  v)j 
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THEORY  OF  THE  STRENGTH  OF  COMPOSITE  PLATES  SYMMETRICALLY  LAMINATFD 

■  *  OFTMOTJPIC'  LAYERS"  AND  SUBJECT'  TO  HF- PLANE'  BIAXIAL  LOADING — 


Unlike  the  theory  of  stiffness  presented  in  Section  2,  there  is  no  universal 
agreement  as  to  the  appropriate  macroscopic  theory  to  describe  the  strength 
of  anisotropic  composites  under  multiaxial  loading.  As  in  the  case  of 
monolithic  metallic  alloys,  certain  theories  fit  certain  materials;  while 
sane  materials  exhihit  different  modes  of  failure  depending  upon  the  temper¬ 
ature,  rate  of  loading,  and  other  environmental  conditions,  and  thus  require 
different  theories  for  different  conditions. 


In  composite  materials,  there  are  at  least  three  distinct  methods  of  pre¬ 
dicting  multiaxial  strength: 

(1)  Empirical  and  semiempirical  macroscopic  theories 

(2)  Fracture  mechanics  (for  brittle  materials) 

(3)  Plastic  tensile  instability  (for  ductile  materials) 

Of  these  three  approaches,  only  the  first  one  is  adequately  developed  at 
present  for  direct  comparison  with  experimental  results.  Approaches  of 
this  type  are  reviewed  in  Section  3a  and  compared  with  the  limited  test 
data  previously  available  in  Section  3b.  The  theories  reviewed  in  Section 
3a  are  applied  to  symmetrically  laminated  composite  plates  in  Section  3c. 


The  application  of  fracture  mechanics  to  anisotropic  materials  is  in  such 
an  embryonic  stage  that  this  approach  has  not  yet  been  extended  to  the  case 
of  multiaxial  loading.  However,  for  those  readers  desiring  information  on 
the  subject  of  fracture  mechanics  of  anisotropic  and  composite  materials, 
this  topic  is  included  in  the  Selected  Bibliography  presented  at  the  end 
of  this  report.  Although  Hill's  theory  of  orthotropic  plasticity  has  been 
applied  to  prediction  of  plastic  tensile  instability  (ductile  fracture)  of 
orthotropic  sheets  subject  to  in-plane  biaxial  loading,  the  Hill  theory 
itself  has  not  been  adequately  verified  for  fibrous  composites.  However, 
theories  of  anisotropic  plasticity  and  plastic  tensile  instability  in  aniso¬ 
tropic  materials  are  included  in  the  Selected  Bibliography. 


a.  Strength  Theories  for  Single-Layer  Qrthotropic  Composite  Materials 

Perhaps  the  most  widely  used  yield- strength  criterion  for  orthotropic 
materials  is  the  following  one  proposed  in  1948  by  Hill  (Reference  34) 
as  follows  (for  the  general  three-dimensional  case): 

H1  *~a2  "  °1-*2  +  H2  ^a3  "  ap2  +  H3  (°i  "  a2^ 2 


+  2H4  o|  +  2  H5  +  2H6  =  1 


(38) 
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Equation  (38)  has  been  found  to  predict  with  sufficient  accuracy  the 
initial  yielding  of  many  metallic  alloys.  Since  Equation  (38)  coincides 
with  the  von  Mises  yield  criterion  for  the  case  of  a  material  which  is 
isotropic,  it  can  be  considered  to  he  a  generalization  of  that  criterion. 
The  major  limitations  of  Equation  (38)  are  that  it  assumes  that  hydro¬ 
static  stress  (Ph  =  °i  +  °2  +  03)  does  not  affect  yielding  and  that  the 
yield  strengths  in  simple  compression  are  identical  to  the  yield  strengths 
in  simple  tension  along  the  same  axes. 

For  the  case  of  a  thin  layer,  a ,  =  o.  =  0,  =  0;  i.e.,  a  state  of  gener¬ 
alized  plane  stress  (GPS)  exists.  Tnen  Equation  (38)  reduces  to  the 
following  relation: 

(°x/ Sx)2  +  (a2/S2)2  -  r’1  (oj/Spoysp  +  (a6/S6)2  =  1  (39) 

It  is  noted  that  Equation  (39)  can  be  visualized  as  a  surface  in  three- 
dimensional  stress  space  (o.,  o2>  o^)  as  shown  in  Figure  2.  Despite  the 
fact  that  the  failure  surface  is  three-dimensional,  the  type  of  failure 
being  depicted  is  due  to  a  planar  stress  system.  This  is  definitely 
more  complicated  than  in  the  case  of  a  simple  isotropic  material,  for 
which  GPS  failure  can  be  represented  by  a  curve  in  two-dimensional  stress 
space  (0.,  o-),  as  discussed  in  References  35  and  36.  The  reason  for 
tiie  additional  complexity  in  the  case  of  materials  which  behave  aniso- 
tropically  is  that  failure  is  definitely  affected  by  the  amount  of  shear 
stress  acting  on  the  planes  normal  to  the  material -symmetry  axes.  In 
contrast,  in  a  simple  isotropic  material,  one  can  select  directions 
(principal-stress  directions)  such  that  no  shear  stress  acts  and  these 
coincide  with  material -symmetry  directions,  since  any  direction  is  one 
of  the  latter.  The  three-dimensional  failure  surface  representation  of 
failure  of  anisotropic  materials  under  GPS  conditions  has  been  used  by 
Norr j  s  (Reference  37) ,  Ashkenazi  (References  38  and  39) ,  and  Chamis 
(Reference  16) .  Finally,  it  should  be  mentioned  that  failure  as  dis¬ 
cussed  in  connection  with  plots  such  as  Figure  3  can  be  defined  in  any 
manner  desired  for  the  particular  application:  initial  yielding,  0.2- 
percent  offset  yield,  100,00  cycles  of  fatigue  loading,  10,000  hours  of 
elevated- temperature  life,  ultimate,  etc.  Here,  failure  is  construed  to 
mean  initial  yielding  only. 

Tsai  (Reference  18)  showed  that  if  the  material  is  assumed  to  be  a  unidi¬ 
rectional  fiber- reinforced  composite  having  the  same  properties  in  the 
thickness  direction  as  it  has  in  the  plane  normal  to  the  fibers,  r  =  S./S2 
and  Equation  (39)  becomes 

(°  1/S1) 2  +  (W2  '  0V2/Sl)  +  (a6/S6)2  =  1  (40) 

Norris  and  McKinnon  (Reference  40)  suggested  the  following  strictly 
empirical  strength  relationship: 

CySj)2  *  (°2/S2)2  +  (°6/S6)2  =  1  (41) 
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Figure  2.  Possible  GPS  Failure  Surface  Predicted  by  Hill  and  Tsai. 
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Later  Norris  (Reference  37)  derived  the  following  relationships  by 
assuming  that  the  distort ional  energy  at  failure  was  a  constant: 

(<y Sj)2  +  Ca2/S2)2  -  (o1/S1)  Ca2/S2)  +  (a^)2  =  1  (42) 

(o2/S2)2  =  1  (43) 

(VSj)2  =  1  (44) 

It  is  noted  that  Equation  (42)  is  a  special  case  (r  =  1)  of  Equation  (39). 
Apparently  Equation  (42)  was  also  derived  independently  by  Ashkenazi 
(Reference  41).  The  GPS  failure  surface  predicted  by  Equations  (42) 
through  (44)  is  shown  in  Figure  3. 

In  several  theories,  provision  was  made  for  hydrostatic-stress  effects; 
namely,  in  the  work  of  Marin  (Reference  42),  Hu  and  Pae  (Reference  43), 
and  Bert  et  al  (Reference  36) .  The  latter  investigators  proposed  the 
following  form  of  relationship  for  the  generalized-plane-stress  case: 

B1  °1  +  B2  °2  +  B3  °1  °2  +  B4  °1  +  B5  °2  =  1  (45) 

It  is  interesting  to  note  that  Equation  (45)  is  the  equation  of  a 
generalized  conic  curve  and  that  it  has  five  constants  which  can  be 
determined  experimentally  by  five  tests  in  the  tens  ion- tens  ion  quadrant: 
uniaxial  tension  at  0  and  90  deg,  and  1:1,  1:2,  and  2:1  biaxial  tension 
at  0  deg.  However,  it  should  be  noted  that  Equation  (45)  is  applicable 
only  to  cases  when  the  principal-stress  directions  coincide  with  the 
material-symmetry  directions,  since  it  does  not  contain  o^. 

The  maximum- strain  failure  criterion  has  been  proposed  by  sane  recent 
investigators  (Reference  44) ;  however,  no  experimental  verification  of 
this  criterion  has  been  published.  Also,  recently  McGill  (Reference  45) 
showed  that  the  classical  representation  of  this  criterion  for  GPS 
conditions  was  incorrect. 

Several  attempts  have  been  made  to  account  for  the  known  difference 
between  tensile  and  compressive  strengths  in  actual  fibrous  composites. 
The  first  of  these  was  Norris,  who  suggested  that  in  applying  Equations 
(42)  through  (44),  if  a.  (i  =  1,  2)  is  tensile,  then  the  tensile  S, 
should  be  used.  Another  approach  was  used  by  Stassi-D’Alia  (Reference 
46) ;  his  proposed  relation  is  as  follows  for  the  case  of  biaxial  princi¬ 
pal  stresses  aligned  with  the  material-synmetry  axes: 

°1  +  °2  ’  °1  °2  +  (K  -  1) (^  +  a2)  =  KsJ  (46) 

C  T  C  T 

where  K  =  S,/S,  and  S,  and  S,  are  the  compressive  and  tensile  strengths 
in  direction  1.  In  tnis  theory,  it  is  assuned  that  yielding  occurs 
when  the  octahedral  shear  stress  reaches  a  value  which  is  proportional 
to  the  hydrostatic  stress.  However,  as  pointed  out  by  Chamis 
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Figure  3.  Possible  GPS  Failure  Surface  Predicted  by  Norris. 

I 


20 


(Reference  16),  in  the  Stassi-D'Alia  theory, it  is  implied  that  the 
following  relationship  holds: 

r  t  f  T 

sys*  =  s~/sj  (47) 


Unfortunately,  actual  filamentary  composites  do  not  generally  satisfy 
Equation  (47) . 

Hoffman  (Reference  47)  has  extended  the  work  of  Hill  to  account  for 
different  strengths  in  tension  and  compression.  For  the  general 
three-dimensional  case,  his  result  can  be  expressed  as  follows: 

H1  (°2  "  °1^2  +  H2  (°3  "  ap2  +  H3  (cti  "  °2 )2  +  2H4  °4 

+  2H5  °5  *  2I,6  “6  *  H7  °1  *  H8  °2  *  «9  °3  '  1  (4« 


The  nine  coefficient  H,  through  Hg  are  uniquely  determined  by  nine 
simple  tests:  three  uniaxial  tensile  strengths,  three  uniaxial  compres¬ 
sive  strengths, and  three  pure  shear  strengths.  When  the  compressive 
strength  along  each  material-symmetry  axis  is  equal  to  the  corresponding 
tensile  strength,  H?  through  Hg  vanish  and  Equation  (48)  reduces  to 
Hill's  relation,  Equation  (38). 


Further  assuming  that  the  material  is  a  unidirectional  one  having  the 
same  properties  in  the  thickness  and  transverse  planar  directions, 
Hoffman  reduced  Equation  (48)  to  the  following  for  the  generalized- 
plane-stress  case: 


(o|/S^)  +  (o2/sJsl[)  -  (o1  o2/sJs[)  +  [(S’JV1 


(49) 

-  (S^)’1]  ox  ♦  [(S^)'1  -  (S^)'1]  a2  +  (a6/S6)2  =  1 
C  T  C  T 

It  is  noted  that  when  and  S~  =  S-,  Equation  (49)  reduces  to 

Tsai's  result,  Equation1 (40) 7  The  GPS  failure  surface  corresponding  to 
Equation  (48)  is  shown  in  Figure  4. 


Gol'denblat  and  Kopnov  (Reference  48)  have  used  tensor  consit  rations  to 
arrive  at  a  failure  theory  even  more  general  than  Hoffman's,  because  they 
took  into  consideration  the  dependence  of  shear  strength  on  direction 
(i.e.,  sign).  ' 


Recently,  Capurso  (Reference  49)  presented  a  yield  condition  which  is  a 
generalization  of  the  Tresca  maximum-shear-stress  yield  condition  to  the 
case  of  an  orthotropic  material  with  different  tensile  and  ccmpressive 
strength.  However,  unlike  all  of  the  other  theories  which  take  into 
consideration  oj  /  0?,  this  theory  is  unaffected  by  hydrostatic  stress. 
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All  of  the  failure  criteria  discussed  so  far  have  been  phenomenological 
in  nature  and  have  not  taken  into  account  any  specific  mechanism  of 
failure  or  any  mechanical  properties  of  the  individual  constituent 
materials.  One  of  the  earliest  attempts  to  take  these  latter  aspects 
into  consideration  was  made  hy  Stowell  and  Liu  (Reference  50) ,  who  pro¬ 
posed  a  simple  mechanistic  model  having  three  failure  modes:  fiber 
failure,  matrix  shear  failure,  and  matrix  transverse  failure.  However, 
they  have  applied  their  criterion  to  only  the  case  of  uniaxial  loading 
applied  at  any  arbitrary  angle  to  the  major  material-symmetry  axis  of  the 
composite. 
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Figure  4.  Possible  GPS  Failure  Surface  Predicted  by  Hoffman's  Criterion 
for  S^/S^  and  Syslj, 

Additional  research  on  the  distort ion- energy  criterion  was  carried  out 
by  Griffith  and  Baldwin  (Reference  51),  and  Chamis  (Reference  16)  extended 
their  work  to  incorporate  different  tensile  and  compressive  strengths. 

For  the  generalized-plane-stress  case,  Chamis'  equation  can  be  written 
as  follows: 
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(°“/S“)2  -  (a*/s|)*  -  Ka@  (c"/|S“|) 
Co|/ | s| | )  ♦  (°6/s6)2  =  1 


(SO) 


where  a,  6  =  T  or  C,  i.e.,  tension  or  compression,  and  K  is  an 
empirical  constant  determined  from  biaxial  tests.  The  wSy  in  whic’ 

Chamis  has  taken  into  consideration  different  tensile  and  compressive 
strengths  is  to  empirically  determine  a  different  K  for  different 
quadrants,  i.e.,  f  K^,  j-  K^.  Also,  as  he  pointed  out,  the  surface 
may  have  discontinuous  derivatives  at  the  quadrant  junctures,  as  shown 
in  Figure  5. 

b.  Comparison  of  Strength  Predictions  with  Existing  Test  Data  for  Single- 
Layer  Orthotropic  Materials 

The  sparse  amount  of  experimental  strength  data  on  composite  materials 
available  for  comparison  with  various  strength  theories  has  been  limited 
to  three  kinds  of  loading:  (1)  uniaxial  loading  (tension  and  compression) 
at  various  off-axis  orientations,  (2)  pure  shear  loading  at  various 
orientations,  and  (3)  biaxial  loading. 

For  the  case  of  uniaxial  loading  at  an  angle  0  with  the  major  material- 
symmetry  axis,  the  applied  uniaxial  stress  a  I  must  be  resolved  into 
components  along  the  material -symmetry  axes  as  follows: 

°1  =  °i  m2 »  °2  =  ai  n2»  a6  =  ”°i  11111  (51) 

where  as  before  m  =  cos  0^  and  h  =  sin  0^. 

These  stress  components  can  now  be  substituted  in  any  desired  strength 
theory  to  obtain  a  theoretical  prediction  of  off-axis  strength. 

The  Forest  Products  Laboratory  has  obtained  fair  to  good  agreement 
between  the  Norris  distortion-energy  strength  theory,  Equation  (42) ,  and 
experimental  uniaxial  data  for  plywood  (Reference  37)  and  various  glass 
fabrics  impregnated  with  various  resins  (References  52,  53,  and  54). 

Recently  Tsai  et  al  (Reference  55)  have  compared  off-axis  uniaxial 
tension  and  compression  data  for  a  unidirectional  glass  fiber/epoxy 
laminate  with  the  predictions  of  the  maximum-normal-stress,  maximum- 
shear-stress,  and  distortion- energy  theories.  They  found  the  best 
agreement  with  the  Tsai  version  of  the  distortion- energy  theory,  Equa¬ 
tion  (40),  where  tensile  properties  (s'*,  Sp  were  used  for  tension 
loading  and  (Sp  sQ)  were  used  for  compressive  loading.  In  Reference 
47,  Hoffman  showed  good  agreement  between  his  generalized  strength  theory 
predictions  and  both  tensile  and  compressive  off-axis  uniaxial  data  for 
GRP. 
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Figure  5.  Possible  GPS  Failure  Surface  Predicted  by  Chamis. 


I 
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For  pure  shear  loading  at  an  off-axis  angle,  the  applied  shear  stress 
a l  must  be  resolved  into  the  following  material- symnetry  stress  compo¬ 
nents  : 

Qj  =  2o£  mn,  =  -2o£  mn,  o^  =  a ^  (m2  -  n2)  (52) 

Data  from  off-axis  pure-shear  tests  have  been  quite  limited.  In  tests 
on  plywood  (Reference  37)  and  various  glass -fabric/resin  combinations 
(References  52,  53,  and  54),  the  agreement  with  the  Norris  distortion- 
energy  criterion  was  rather  poor.  Only  very  limited  data  are  available 
on  strength  of  composites  under  biaxial  loading,  apparently  all  for  the 
aligned  case  (i.e.,  0=0).  Ely  (Reference  56)  tested  glass -fabric/epoxy 
circular  tubes  subject  to  combinations  of  hoop  tension  and  axial  compres¬ 
sion.  He  used  the  Stassi-D'Alia  theory  to  test  his  results,  with  good 
agreement.  However,  Chamis  (Reference  16)  showed  later  that  Ely's  test 
results  could  also  be  represented  by  the  Chamis  empirical  criterion. 

Weng  (Reference  57)  tested  two  grades  of  graphite  composite, in  the  form 
of  circular  tubes,  under  biaxial  loading  in  all  four  quadrants  (tension 
T/compression  C,  T/T,  C/T,  C/C).  The  results  were  fitted  well  bv  the 
Chamis  empirical  criterion  with  a  different  value  of  K  for  eacn  quad¬ 
rant  (see  Table  II .)  P 


TABLE  II.  VALUES  OF 
COMPOSITE 

A$) 

IN  THE  CHAMIS  EMPIRICAL  CRITERION  FOR  VARIOUS 
BIAXIAL- LOADING  QUADRANTS 

Material 

Quadrant 

Glass -fabric/ epoxy 

C/T 

-  0.36 

Grade  JT  graphite 

T/T 

0.50  to  0.90 

T/C 

0.40 

C/C 

l.QQ  to  1.45 

C/T 

0.40 

Grade  ATJ  graphite 

T/T 

-  0.6  to  -  1.60 

T/C 

0  to  1.20 

C/C 

0  to  -  1.20 

C/T 

0.6  to  1.20 
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Gol'denblat  and  Kopnov  (Reference  48)  showed  good  agreement  with  their 
theory  for  two  kinds  of  FRP  in  all  four  quadrants. 

In  the  present  work,  same  of  the  biaxial- load  tests  were  conducted  with 
the  biaxial  principal  stresses  a'  and  o ^  aligned  at  an  angle  (e  -  45  deg) 
to  the  material- symmetry  axes,  in  this  case,  the  stress  components 
associated  with  the  material-symmetry  axes  can  be  calculated  from  the 
following  equations: 


0  =  oj  m2  +  a*  h2 

°2  =  +  °2  m2  (53) 

°6  =  '  (ai  "  ^ 

c*  Strength  of  Symmetrically  Laminated  Composite  Plates 

The  strength  of  laminated  composite  plates  can  be  quite  different  than 
that  of  a  single  layer,  even  in  the  case  of  symmetrical  lamination. 

This  is  due  to  three  factors:  (1)  the  heterogeneity  of  the  composite 
in  the  thickness  direction,  (2)  the  thermally  induced  residual  stresses, 
and  (3)  interply  strain  energy. 

Perhaps  the  first  investigator  to  consider  the  heterogeneity  of  the 
composite  in  the  thickness  direction  was  Tsai  (Reference  18) ,  who  applied 
a  failure  criterion  to  each  ply  to  determine  which  one  failed  (yielded) 
first  and  then  used  degraded  properties  for  each  failed  ply  to  determine 
a  new  stress  distribution  until  complete  failure  (i.e.,  failure  of  all 
of  the  plies)  occurred.  This  concept  was  extended  further  by  Tsai  et  al 
(Reference  55) . 

Tsai  (Reference  18)  also  included  the  effect  of  the  residual  stresses 
which  are  induced  thermally  due  to  the  difference  between  the  laminate 
curing  temperature  and  the  operating  temperature.  Since  this  factor  is 
very  nearly  constant  in  the  experimental  program  reported  here,  as  men¬ 
tioned  in  Section  2b,  this  effect  is  emitted  in  the  present  analytical 
work. 

The  effect  of  interply  shear  strain  energy  on  strength  was  discussed  in 
detail  by  Chamis  (Reference  16).  However,  as  mentioned  in  Section  2b, 
this  effect  is  relatively  small  and  thus  is  neglected  in  the  present 
analysis. 

Although,  so  far  as  is  knewn  by  the  present  investigators,  the  effect  of 
filament  crossovers  has  not  been  considered  theoretically  in  any  pub¬ 
lished  analysis,  it  was  discussed  qualitatively  by  Tsai  et  at  (Reference 
55) . 

For  the  present  case,  the  use  of  symmetric  lamination  arrangements 
eliminates  coupling  between  in-plane  and  out-of -plane  effects. 
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Furthermore,  since  the  loading  is  applied  in  the  plane,  no  bending 
moments  are  induced.  Under  these  simplified  conditions,  Tsai’s  general 
analysis  is  reduced  considerably.  The  strains  can  he  found  by  invert¬ 
ing  Equation  (19)  as  follows: 


XX 

A11 

A12 

A16 

N 

XX 

yy 

8 

A12 

A22 

A26 

N 

yy 

(54) 

*y_ 

_A16 

A26 

A66_ 

>J 

where  the  [A*.]  matrix  represents  the  inverse  of  the  [A-.]  matrix.  An 
alternative  out  equivalent  approach  is  to  use  Equations1^ 23)  to  compute 
Eh,  ^22’  v12  t0  use  f°ll0WinS  relationships: 


exx 

(%)  1  _vi2/Fii 

Nxx 

Eyy 

=  (l/h) 

■v12/E11  ^22*  ^2p 

N 

yy 

exy 

Cey'1  (ey-1 

LxyJ 

where 


(H'il)'1  =  (Fn)"1  m4  +  L  m2  n2  +  (E^)'1  n4 

(F’2 )_1  =  (Fn)’1  n4  +  r  m2  n2  +  (E^)’1  m4 

^12/Fil  =  (°12/E lP  -  ®  ™ 

(C' ) "  *  =  (G)  ’  *  +  4£T  mn 

(ei2)_1  =  [(E22)_1  n2  -  (En)'1  m2  +  (L/2)(m2  -  n2)]  2  mh 
(f ijj)"1  =  [(E22)_1  m2  -  (En)'1  n2  +  (L/2) (m2  -  n2)]  2  mh 
L  =  (G)  -  (vi2/En)  -  (v2]/^22^ 

FI  =  (En)_1  +  (F22)_1  -  r 


Once  the  strains  are  determined  from  either  Equation  (54)  or  Equation 
(55) ,  the  stresses  in  the  respective  layer  can  be  calculated  by  use  of 
Equation  (16) .  Then  the  stress  components  can  be  expressed  in  terms 
of  the  desired  strength  theory  as  presented  in  Section  3a. 
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SPECIMEN  DESIGN,  DEVELOPMENT, 


AND  PRELIMINARY  EVALUATION 


As  mentioned  in  Section  1,  the  primary  objective  of  this  research  is  the 
mechanical  characterization  of  laminated. composite  materials,  with  parti¬ 
cular  emphasis,  on  biaxial  loading.  The  separate  determination  of  the 
uniaxial  tensile  and  shear  properties  was  necessary  to  complete  the  charac¬ 
terization  of  the  material  and  to  permit  correlation  with  theoretical  pre¬ 
dictions  of  stiffness  and  strength. 

This  section  discusses  the  basic  philosophy,  experimental  development,  and 
final  design  of  the  test  specimens  in  the  data  gathering  portion  of  the 
research  program. 

a.  Biaxial-Loading  Specimens 

Reference  36  contains  a  critique  of  the  various  kinds  of  specimens  used 
to  obtain  biaxial-load  test  data  for  metallic  alloys.  Some  of  the  same 
kinds  of  specimens  have  been  used  to  obtain  limited  biaxial- load  test 
data  for  composites.  Closed-end,  thin-walled,  internally  pressurized 
cylinders  of  glass-filament -wound  construction  have  been  used  to  simu¬ 
late  solid-propellant  rocket-motor  casings  (Reference  58).  However, 
only  one  biaxial-load  condition  (tension  1:2)  was  covered,  and  the 
results  of  such  tests  are  affected  by  stress  nonuniformities  near  the 
ends.  To  obtain  data  for  deep-diving  submersible  applications,  sane 
tests  have  been  conducted  on  cylindrical  shells  subject  to  external 
pressure,  giving  a  1:1  canpression  biaxial-load  condition  (Reference  59). 

Recently  cruciform  sandwich-beam  specimens  have  been  used  to  deteimine 
biaxial  yield  strengths  of  boron- fiber- reinforced  plastic  (Reference  60). 
The  disadvantages  of  this  type  of  specimen  are  the  stress  concentrations 
produced  at  the  comer  fillets,  the  possible  effect  of  the  core-to- 
facing  bond  on  the  laminate  strength,  nnu  the  fact  that  the  facings  expe¬ 
rience  some  bending  action  in  addition  to  the  primary  membrane  action. 

In  biaxial-load  tests  on  metals,  Terry  and  McClaren  (Reference  61)  used 
flat  specimens  with  a  reduced- thickness  test  section  and  four  loading 
pins.  In  crack  propagation  tests  under  biaxial  loading,  Douglas  Aircraft 
Division  (Reference  62)  used  a  machined  slit  to  simulate  a  crack  (instead 
of  a  reduced  thickness  test  section)  and  eight  loading  pins.  In  sane 
very  recent  FRP  laminates  tests  reported  after  the  present  program  began, 
there  was  no  reduced- thickness  section  and  four  loading  pins  were  used 
(Reference  63) . 

The  basic  philosophy  underlying  the  development  of  the  biaxial- loading 
specimens  was  to  evolve  flat  laminate  specimens  which  were  simple  to 
fabricate  and  yet  capable  of  developing  high- stress  levels  and  a  uniform 
tension  stress  field  of  prescribed  biaxial- load  ratio  in  the  test 
section. 

To  eliminate  bending  and  twisting  effects  induced  by  in-plane  loading, 
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it  was  decided  to  use  symmetrical  laminating  arrangements.  This  eliminated 
coupling  between  in-plane  and  out-of-plane  effects,  as  discussed  in  Section 


Preliminary  tests  indicated  that  the  thicker  section  of  the  specimens 
should  be  approximately  three  times  as  thick  as  the  reduced- thickness 
test  section.  This  choice  is  predicated  on  the  following  considerations: 

If  there  is  not  sufficient  difference  in  thickness  between  the  test  sec¬ 
tion  and  the  remainder  of  the  specimen,  failure  would  usually  occur  in 
the  vicinity  of  the  grips.  (This  occurred  in  all  of  the  nonreduced- 
thickness- section  specimens  reported  in  Reference  63.)  On  the  other  hand, 
if  there  is  too  much  difference  in  thickness,  three-dimensional  geome¬ 
tric  stress  concentration  would  artificially  raise  the  stress  level  at 
the  edge  of  the  test  section.  Fortunately,  the  natural  filleting  which 
occurs  in  the  fabrication  of  the  specimens  (see  Section  6a)  helps  to 
reduce  this  stress  concentration.  Another  limitation  is  the  maximum 
thickness  of  laminate  from  which  excess  resin  can  be  successfully  bled. 

A  series  of  preliminary  uniaxial  tension  tests  brought  out  what  turned 
out  to  be  the  most  critical  problem  in  the  development  of  the  specimens: 
the  problem  of  achieving  adequate  load- carrying  capacity  in  the  metal 
tabs  used  to  transfer  the  load  fran  the  loading  device  to  the  specimen. 

(The  loading  device  is  described  in  Section  5a.)  Interspersing  layers  of 
m^allic  shim  stock  between  the  FRP  layers,  clanping  the  tabs  exter¬ 
nally,  and  adding  a  ceramic  thickener  to  the  bonding  adhesive  did  not 
help.  Even  with  the  use  of  tapered- thickness  tabs,  a  shear  load-carrying 
capacity  of  only  1500  psi  was  achieved.  For  simplicity,  a  minimum  number 
of  loading  pins  (four)  was  used  in  the  biaxial- loading  specimens. 

In-plane  biaxial- load  tests  have  previously  been  carried  out  on  metals 
by  Terry  and  McClaren  (Reference  61) .  They  stated  that  after  trying 
various  designs  for  the  contour  of  the  reduced- thickness  test  section, 
they  selected  their  final  design,  which  was  a  square  contour  with 
rounded  corners.  Although  it  was  stated  that  uniformity  of  strain  distri¬ 
bution  was  the  criterion  used  for  selection  of  the  final  contour  design, 
no  strain  distribution  data  were  reported. 

The  test-section  contours  used  in  the  present  program  were  based  on  two- 
dimensional  elasticity  theory,  as  descrihed  in  Appendix  I .  The  resulting 
contours  were  a  circular  contour  for  the  specimens  subject  to  a  biaxial- 
load  ratio  of  1:1  and  a  1:/?  elliptical  contour  for  the  hiaxial-load  ratio 
of  1:2.  It  should  be  noted  that  the  reduced- thickness  contour  synthesized 
in  Appendix  I  results  theoretically  in  a  uniform  biaxial-stress  state  in 
the  test  section  only.  To  achieve  a  uniform  stress  state  in  the  surrounding 
thicker  section  would  require  a  local  ring- type  reinforcement  of  varying 
cross-sectional  area.  However,  it  is  more  desirable  to  have  a  nonunifonn 
stress  field  in  the  thicker  section:  higher  stresses  adjacent  to  the  test 
section,  gradually  decreasing  as  the  distance  from  test  section  is  increased 
(to  compensate  for  the  increased  localized  stresses  in  the  vicinity  of  the 
loading  tabs) . 
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The  scale  (i.e.,  the  actual  size)  of  the  test-section  contour  was 
based  on  the  following  considerations:  The  test  section  must  he  large 
enough  to  allow  an  ample  sample  of  material  and  in  particular  to  he 
sufficiently  larger  than  the  strain  gages  used.  On  the  other  hand, 
space  limitations  on  the  maximum  size  of  the  loading  fixture  and  provi¬ 
sion  for  adequate  distance  between  the  edges  of  the  specimen  and  the 
edge  of  the  test  section  placed  a  maximum  limit  on  the  test-section 
size. 

Preliminary  tests  on  1:1  biaxial  specimens  indicated  that  a  minimum 
test-section  thickness  of  four  plies  was  necessary  to  minimize  specimen 
warpage.  Also,  in  these'  preliminary  tests,  brittle  lacquer  was  applied 
to  the  specimen  to  obtain  a  qualitative  indication  of  the  stress  field 
throughout  the  specimen.*  Over  the  whole  reduced- thickness  test  section, 
the  brittle- lacquer  crack  pattern  was  bidirectional,  like  that  showr  in 
Reference  64;  this  indicated  that  a  biaxial  stress  field  was  developed 
throughout  the  test  section.  In  the  thicker  section  adjacent  to  the 
test  section,  the  cracks  were  radially  oriented,  thus  indicating  that 
the  maximum  stresses  were  circumferential.  These  brittle- lacquer  tests 
indicated  a  very  highly  stressed  region  at  the  edge  of  the  specimen 
between  two  adjacent  loading  tabs.  To  alleviate  this  undesirable  condi¬ 
tion,  the  sides  of  loading  tabs  were  rounded,  and  considerably  more 
care  was  taken  in  achieving  smooth,  generously  curved  outer  edges  of 
the  specimens. 

Additional  preliminary  strain  evaluations  were  carried  out  on  a  1:1 
biaxial  specimen  using  metallic-foil- type  electric- resistance  strain 
gages  (Budd  C6-141B  single  gages  and  Budd  C6-141-R3V  strain  rosettes). 

The  gage  locations  are  shown  in  Figure  6,  and  the  load-strain  data  are 
presented  in  Figure  7.  The  good  uniformity  of  the  strain-gage  readings 
associated  with  the  same  orientation  in  the  test  section  substantiates 
the  validity  of  the  specimen  design.  The  slight  difference  due  to 
orientation  is  attributed  to  the  slight  difference  between  E,,  (in  the 
warp  direction)  and  E22  (in  the  weave  direction)  for  this  material.  As 
a  result  of  the  straifTunifoimity  with  respect  to  orientation  within 
the  test  section,  it  was  decided  to  use  only  two  strain  gages  in  the 
main  data- gathering  tests  rather  than  a  more  expensive  three- element 
rosette  which  would  be  necessary  if  neither  the  principal-stress 
direction  nor  the  principal- strain  direction  were  known. 

Since  the  highest  strains  were  recorded  in  the  region  between  two 
adjacent  loading  tabs,  the  loading- tab  external- comer  radii  were 
increased  to  the  final  value  of  1.5  inches.  The  final  proof  of  the 


*Although  photoelastic  coatings  used  in  conjunction  with  a  reflective 
polariscope  are  generally  more  accurate  than  brittle  lacquers,  it  was  not 
possible  to  find  a  photoelastic  coating  having  a  satisfactory  combination 
of  high  stress-optical  sensitivity  and  low  in-plane  stiffness  (product  of 
Young's  mudulus  and  coating  thickness). 
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STRAIN  -  INCHES /TNCH  X  10' 
a.  Gages  1  Through  9 

Load-Strain  Data  on  Preliminary  Strain  Field  Evaluation  of  1:1 
Biaxial-Loading  Specimen. 


STRAIN  -  INCHES/INCH  X  10 
b.  Gages  10  through  21 


Figure  7.  Continued 


validity  of  the  tests  in  this  regard  was  that  the  fracture  went  through 
at  least  a  part  of  the  test  section  in  36  of  the  54  final  data- gathering 
specimens.  (This  compares  very  favorably  with  the  zero  percent  failures 
in  the  main  section  in  the  tests  reported  by  Grimes  et  al  in  Reference 
63) . 

In  view  of  the  good  strain  distribution  achieved  in  the  1:1  biaxial- 
load  specimens,  a  much  more  limited  strain- field  survey  was  conducted 
for  the  1:2  biaxial-load  specimens  (see  Figure  8). 

The  final  specimen  designs  are  shown  in  Figure  9.  The  strain  gages 
used  were  metallic-foil  type  (Budd  C6-141B),  mounted  as  shown  for  the 
main  data-gathering  tests. 

b.  Torsion-Tube  Shear  Specimens 

There  has  been  considerable  controversy  regarding  the  best  way  to  obtain 
pure  shear  data  for  composite  materials.  The  in-plane  panel  shear  or 
picture- frame  method  has  been  used  most  extensively  in  the  past  (Refer¬ 
ences  52  and  65).  However,  this  method  has  been  subject  to  criticism 
because  the  stress  field  deviates  substantially  from  pure  shear  (Refer¬ 
ence  66) . 

The  off-axis  tensile  test  has  been  used  in  conjunction  with  orthotropic 
elasticity  relationships  to  determine  the  shear  modulus  (Reference  54). 
However,  this  method  has  been  questioned  recently,  since  off-axis  load¬ 
ing  of  anisotropic  materials  induces  in-plane  bending  (Reference  67)  and 
since,  due  to  the  cut-fiber  effect,  the  test  is  thought  to  be  more  of  an 
in-plane  shear  test  of  the  resin  only  (Reference  9).  The  so-called  bow- 
tie  specimen  was  devised  to  eliminate  both  of  these  objections;  however, 
it  also  has  not  been  entirely  satisfactory  due  to  the  stress  concentra¬ 
tion  in  the  fillet  at  the  minimum  cross  section  (Reference  68) . 

Some  of  the  tests  used  to  determine  shear  properties  of  composites 
result  in  out-of-plane  shear  loading,  i.e.,  out-of-plane  twisting,  and 
thus  are  not  valid  for  determining  the  in-plane  shear  modulus.  One 
widely  used  test  of  this  type  is  the  so-called  plate  twist  test  in  which 
a  square  plate  is  loaded  downward  at  two  diagonal  comers  and  upward 
at  the  other  two  diagonal  comers  (Reference  69) .  Another  is  the  split¬ 
ring  out-of-plane  shear  test  (Reference  70) . 

Many  other  types  of  shear  specimens,  such  as  double-notched  and  S-shaped 
tension  specimens  (Reference  71),  deep,  short-span  beams  (Reference  72), 
and  solid  circular  torsion  bars  (Reference  73) ,  have  been  investigated 
to  a  limited  extent.  However,  the  most  popular  shear  test  at  present  is 
the  torsion  of  a  thin-walled  cylinder.  Whitney  (Reference  74)  has  shown 
that  such  a  test  is  insufficient  to  obtain  the  shear  modulus  for  unsym- 
metrically  laminated  composites.  However,  in  the  present  research,  only 
symmetrical  lamination  arrangements  are  used;  thus,  the  test  is  suffi¬ 
cient  to  determine  the  laminate  shear  modulus  directly.  The  tors  ion- tube 
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9.12  Typical  (Horiz . &Vert . 


All  dimensions  Are  in  Inches 


Figure  9.  Biaxial -Loading  Specimen  Designs. 


shear  test  is  believed  to  be  the  only  method  which  meets  all  of  the 
following  requirements  for  a  good  shear  test  method: 

(1)  Development  of  a  uniform  state  of  pure  shear  stress  over 
a  significant  portion  of  the  specimen. 

(2)  Capability  of  determining  both  shear  modulus  and  shear 
strength  in  one  test  of  a  single  specimen  with  a  minimum 
of  data  reduction. 

(3)  Ease  of  specimen  fabrication  from  a  minimum  amount  of 
pre impregnated  laminate. 

(4)  Simplicity  of  test  apparatus. 

(5)  Simplicity  and  repeatability  of  test  procedure. 


Preliminary  torsion  tests  were  conducted  on  tubes  laminated  of  two  plies 
oriented  longitudinally,  to  check  out  the  torsion  test  equipment  (see 
Section  5b).  The  dimensional  and  failure  data  for  the  specimens  are 
listed  in  Table  III. 


TABLE 

III.  RESULTS  OF 

PRELIMINARY  TORSION-TUBE  TESTS_ 1 

Nominal  Diameter, 
in. 

Operating 
Length,  in. 

Shear  Stress  at 
Failure,  psi 

Type  of 
Failure 

Calculated 

‘V’cr’  Psi 

1.590 

1.907 

18,100 

Torsional 

buckling 

18,600 

1.590 

1.250 

28,000 

Torsional 

buckling 

23,000 

Since  these  two  specimens  failed  due  to  torsional  buckling,  calculations 
were  made  using  the  following  approximate  equation,  derived  from  the 
work  of  Simitses  (Reference  75),  for  the  critical  shear  stress 
for  torsional  buckling  of  orthotropic  tubes  with  clamped  edges:  ™ 


(a 


xy-^cr 


(1.03  tt2/12) 


E  5/8 
c22 


EjJ/8  tS/4  l-1/2  r3/4 


(56) 


where  t  =  wall  thickness,  L  =  operating  length,  and  R  =  mean  radius.  As 
can  be  seen  in  Table  III,  the  agreements  between  the  buckling  stresses 
and  the  theoretical  predictions  were  fairly  good. 


To  ensure  that  the  undesirable  buckling  would  not  occur  in  the  final 
specimens,  on  the  basis  of  the  above  equation,  it  was  decided  to  use 
an  operating  length  of  only  5/8  inch  and  an  inside  diameter  of  1.50 
inches. 
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The  use  of  such  a  short  operating  length,  raised  several  questions.  One 
was  the  possibility  of  end  effects  on  the  shear- stress  distribution. 

This  was  investigated  theoretically,  and  it  was  found  that  a  thin-walled 
tuhe  clamped  at  the  ends,  but  free  to  move  axially  and  subject  to  torsion 
loading,  can  experience  only  axial  and  circumferential  displacement  (no 
normal  displacement)  so  long  as  the  huckling  load  is  not  reached. 

Another  question  was  the  possibility  of  an  effect  of  length  on  the  shear 
strength  of  the  composite,  since  increasing  the  length  of  single  glass 
fibers  decreases  their  tensile  strength  (Reference  76).  This  question 
was  investigated  by  testing  three  tubes,  each  having  parallel  plies 
oriented  longitudinally. 

One  specimen  had  an  operating  length  of  3  inches  and  six  plies,  the 
second  had  a  length  of  5/8  inch  and  six  plies,  and  the  third  had  a 
length  of  5/8  inch  and  two  plies.  All  of  the  specimens  had  very  nearly 
the  same  shear  stress-strain  results,  and  thus  it  was  decided  to  stan¬ 
dardize  on  the  operating  length  of  5/8  inch  for  all  of  the  torsion 
tubes.  The  final  specimen  dimensions  and  strain-gage  locations  are 
shown  in  Figure  10. 

c.  Uniaxial  Tension  Specimens 

In  previous  research  on  FRP  laminate  properties  (Reference  1),  the 
tensile-test  specimen  design  selected  was  a  straight-sided  specimen 
with  very  shallow  circular  notches  (0.007  inch  deep)  on  each  side  and 
having  a  net-section  width  of  0.75  inch.  The  purposes  of  the  notch 
were  to  minimize  failures  in  the  testing-machine  grips  and  to  cause 
failure  to  occur  within  the  gage  length  of  the  extensometer,  without 
significantly  reducing  the  strength  values  due  to  geometric  stress 
concentration.  Since  the  notch  was  very  mild  and  the  extensometer 
gage  length  was  appreciable,  the  notch  had  a  negligible  effect  on  the 
extensometer  data.  However,  the  notch  undoubtedly  had  seme  effect  on 
reducing  the  ultimate  strength. 

In  the  present  research,  in  addition  to  strength  and  axial  strain  data 
(to  calculate  the  longitudinal  Young's  modulus  E^) ,  transverse  strain 
data  were  required  (to  calculate  Poisson's  ratio  v,2) .  Since  it  was 
concluded  by  Trantina  (Reference  77)  that  strain  gages  can  be  used  to 
obtain  reliable  elastic-coefficient  data  for  fiher- reinforced  composites, 
preliminary  tensile  tests  were  conducted  on  parallel-ply  specimens 
identical  to  those  described  in  Reference  1,  except  that  a  longitudinal 
and  a  transverse  strain  gage  were  mounted  at  the  center  of  the  specimen. 
The  results  were  unsatisfactory  because  the  notches  increased  the  longi¬ 
tudinal  strain  proportionally  more  than  they  did  the  transverse  strain. 

In  fact,  the  data  were  not  in  agreement  with  Equation  (10),  which  is  a 
fundamental  relation  based  on  basic  thermodynamic  considerations  (see 
Section  2a). 

To  alleviate  the  difficulty  described  above,  plain,  straight- sided 
specimens  (see  Figure  10),  instrumented  with  two  strain  gages,  were 
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Small  semicircular  notch  (0.09  inch  radius) 
made  after  strain-gage  data  were  taken 


b.  Uniaxial  Tension  Specimens 

Figure  10.  Shear  and  Uniaxial  Tension  Specimens. 
(All  Dimensions  Are  in  Inches.) 
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used  to  obtain  the  strain  data.  Then  the  specimens  were  notched  with 
shallow  circular  notches  (0.Q1  inch  deep)  and  tested  to  failure. 

The  question  arose  as  to  whether  the  biaxial  strength  values  differed 
from  the  uniaxial  ones  due  to  the  differences  in  specimen  thickness, 
length,  and  width  between  the  two  types  of  specimens.  To  help  provide 
an  answer  to  this  question,  two  additional  series  of  specimens  were 
made  and  tested.  One  series  of  specimens  was  designed  to  investigate 
the  effects  of  number  of  plies  and  the  length  between  grips,  with  speci¬ 
men  width  held  constant  at  the  same  value  as  the  specimens  shown  in 
Figure  10.  The  test  results  for  the  parallel- laminated  specimens  loaded 
at  0  deg  are  shown  in  Figure  11,  while  Figure  12  gives  limited  test 
results  for  the  parallel-ply  90-deg  specimens,  cross-ply  specimens,  and 
quasi -isotropic  specimens.  The  decrease  in  strength  with  increasing 
test-section  length  follows  the  trend  reported  by  Metcalfe  and  Schmitz 
(Reference  41)  for  glass  fibers  alone,  and  the  increase  in  strength 
with  increasing  nunber  of  plies  is  in  agreement  with  FRP  test  results 
reported  by  many  laboratories,  and  recently  explained  on  a  quantitative 
statistical  basis  by  Scop  and  Argon  (Reference  78). 

The  second  series  of  additional  specimens  was  to  study  effects  of  speci¬ 
men  width  and  stepped  thickness.  Seme  specimens  were  simple  rectangular 
but  with  a  length  equal  to  the  distance  between  the  loading  tabs  on  the 
biaxial  specimens.  As  was  expected  on  the  basis  of  width-effect  tests 
reported  by  Kaman  Aircraft  Corporation  (Reference  79) ,  the  strength 
decreased  with  increasing  width  (see  the  circle  data  points  in  Figure 
13).  The  remainder  of  the  specimens  were  cut  from  some  1:1  biaxial 
specimens  to  various  widths  and  loaded  uniaxially  with  9.25  inches 
between  grips.  Thus,  these  specimens  had  a  reduced  thickness  with  a 
step  change  in  thickness.  The  results  of  these  tests  are  shown  in 
Figure  13  as  square  data  points.  Surprisingly  these  specimens  exhibited 
an  increase  in  strength  with  increasing  width.  This  opposite  effect  is 
obviously  due  to  the  different  restraint  in  the  stepped- thickness 
specimens . 
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ULTIMATE  TENSILE  STRENGTH  -  PS I 


ULTIMATE  TENSILE  STRENGTH  -  PSI  X  10 
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Figure  12.  Effect  of  Number  of  Plies  and  Length  Between  Grips  on  the 

Strength  of  Rectangular,  Uniform-Thickness  Spec mens,  0.750 
Inch  Wide,  With  Various  Lamination  Arrangements  and  Loading 
Directions. 
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ULTIMATE  STRENGTH  -  PSI  X  10 


5 .  SPECIAL  TEST  EQUIPMENT 

The  design  and  operational  features  of  the  special  test  equipment  developed 

in  this  research  program  are  described  in  this  section. 

a.  Biaxial-Loading  Fixture 

In  previous  in-plane  biaxial- loading  fixtures,  two  or  more  hydraulic 
rams  were  used  to  load  the  specimen  (References  61  and  62) .  However, 
it  is  a  difficult  problem  in  hydraulic  control  to  maintain  the  same 
biaxial- load  ratio  throughout  the  test  of  an  anisotropic  composite.  To 
eliminate  this  control  problem,  in  the  present  program,  an  original 
cable- and-pulley  system  reacted  by  a  rigid  steel  frame  was  used  to  load 
the  specimen  biaxially.  Simplicity  is  the  major  feature  of  this  loading 
device. 

After  the  present  device  was  developed,  the  authors  learned  of  the  devel¬ 
opment  of  another  simple  biaxial- loading  device  by  Grimes  et  al  (Refer¬ 
ence  63).  Although  their  device  is  probably  simpler,  it  developed  only 
one  biaxial-load  ratio  and  their  test  results  were  somewhat  questionable, 
as  discussed  later. 

Figure  14  is  a  schematic  diagram  of  the  basic  concept  of  the  cable-and- 
piulley  biaxial- loading  fixture.  It  is  noted  that  the  fixture  can  be 
converted  from  a  1:1  biaxial- load  ratio  to  a  1:2  ratio  by  merely  changing 
the  cable  path  and  pulleys  used.  Figure  15  shows  the  fixture  in  place 
between  the  platens  of  the  60, 000-pound- capacity  Riehic  hydraulically 
operated  universal  testing  machine. 

Although  the  pulley  system  used  in  the  biaxial- loading  fixture  appears 
at  first  glance  to  be  unduly  complicated,  it  was  gradually  evolved  by 
means  of  a  wooden  mockup  and  found  to  be  most  advantageous  in  meeting 
these  design  objectives:  (1)  development  of  maximum  load  in  the  specimen 
equal  to  twice  that  applied  by  the  testing  machine  in  the  case  of  both 
1:1  and  1:2  loading,  (2)  minimization  of  the  maximum  load  applied  to  any 
bearing  in  the  fixture,  (3)  sufficiently  compact  overall  dimensions  to 
fit  between  the  platens  of  the  universal  testing  machine,  (4)  ability 
to  change  biaxial- load  ratio,  and  (5)  allowance  for  adequate  space  for 
the  specimen  and  loading  tabs. 

The  cable  type  and  size  were  chosen  to  meet  these  two  objectives:  (1) 
maximum  flexibility  to  enable  easy  bending  over  the  pulleys  during  opera¬ 
tion,  and  (2)  adequate  static  load  capacity.  Since  the  motion  involved 
during  loading  was  very  small,  the  cable  was  sized  on  the  basis  of  a 
strictly  static  load  (no  motion  over  the  pulleys).  The  cable  finally 
selected  was  5/'8-inch  6  x  36. 

The  standard  cable-socket  method  was  used  to  attach  to  the  cable  for 
maximum  load  capacity.  This  method  entailed  putting  the  cable  end  into 
a  tapered  steel  socket  and  flaring  the  cable  end  until  each  strand  was 
separated.  Next  the  individual  strand  ends  were  bent  back  approximately 
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Figure  15.  Biaxial  Device  Installed  in  Universal  Test  Machine. 
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90  to  180  deg.  Then  a  high- zinc  No.  1  Babbitt  was  poured  into  the  socket, 
forming  a  solid  mass  which  resulted  in  a  high  load-carrying  capability  as 
verified  by  preliminary,  tensile  tests  before  cable  was  installed  in  the 
loading  fixture. 

The  pulley  size  was  standardized  at  a  diameter  of  6  inches  for  all  of 
the  pulleys.  This  size  was  much  smaller  than  that  recommended  by  the 
cable  manufacturer;  however,  it  proved  to  be  adequate  in  this  negligible- 
movement  application.  Each  pulley  contained  a  cylindrical -roller  bearing, 
which  operated  on  a  1.5 -in. -diameter  hardened  steel  shaft,  which,  in  turn, 
was  pressed  into  the  loading  fixture. 

To  permit  the  specimen  to  be  self-aligning,  a  load  link  was  provided 
between  each  loading  tab  on  the  specimen  and  the  biaxial -loading  fixture. 

The  loading  tabs  by  which  the  biaxial  loads  were  transferred  from  the 
load  link  in  the  biaxial -loading  fixture  to  the  specimen  were  made  of 
hardened  SAE  4140  steel.  Eight  tabs  were  used  for  each  specimen:  one 
on  each  side  of  the  specimen  (to  minimize  out -of -plane  bending  in  the 
specimen)  at  each  of  the  four  loading  ears.  To  provide  the  necessary 
shear  area,  the  tabs  were  wider  (7.81  inches  wide)  where  they  were  bonded 
to  the  specimen.  To  promote  a  more  uniform  transfer  of  load  from  the 
steel  loading  tabs  to  the  FRP  specimen,  the  thickness  of  the  tabs  was 
tapered  to  a  minimum  of  0.0625  inch. 

Prior  to  testing  FRP  specimens  in  a  given  configuration,  the  cables  of 
the  loading  fixture  were  "set"  by  applying  load  to  an  aluminum  plate. 

b .  Torsion-Tube  Shear  Test  Apparatus 

To  permit  proper  testing  of  the  tubular  specimens  (see  Section  4b) , 
rather  massive  end  fittings  were  designed  (see  Figure  16).  One  of 
these  was  bonded  to  each  end  of  the  specimen  over  a  length  of  1.5  inches 
in  order  to  achieve  adequate  shear  load-carrying  capacity.  The  end 
fittings  were  mounted  on  a  shaft,  and  one  fitting  was  free  to  slide 
axially  on  the  shaft  to  prevent  the  development  of  axial  forces  as  a 
consequence  of  change  in  specimen  length  during  loading. 

Pure  torque  was  applied  to  the  specimen  via  the  end  fittings  by  means  of 
a  protrusion  on  each  end  fitting.  The  end  fittings  were  bonded  to  the 
specimen  in  such  a  way  that  the  protrusions  were  180  degrees  apart.  This 
minimized  the  change  in  moment  arm  as  the  specimen  twisted  during  loading. 
However,  the  total  angular  rotation  to  failure  was  very  small  in  all  cases, 
since  the  operating  length  of  the  tube  was  quite  short. 

The  loading  protrusions  on  the  end  fittings  were  loaded  by  another  fitting 
on  the  end  of  the  hydraulic  ram.  Figure  17  is  a  photograph  of  the  shear 
test  apparatus  with  specimen  bonded  in  place. 
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Figure  17.  Shear  Test  Apparatus  With  Specimen  Bonded  in  Place. 
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6.  EXPERIMENTAL  PROCEDURES 


a.  Specimen  Fabrication 

All  of  the  specimens  were  fabricated  by  OURI  from  FRP  laminates  consist¬ 
ing  of  181- style-weave,  E-glass  cloth,  with,  a  Volan-A-finish  and  impreg¬ 
nated  with  Epon-828  epoxy  activated  by  curing  agent  Z.  First  the 
material  wus  made  in  the  form  of  a  two-  or  three- layer  preimpregnated 
laminate  (called  prepreg),  using  the  wide  multilayer  prepreg  machine 
reported  in  Reference  9.  Prepreg  52  inches  wide  and  in  lengths  up  to 
2800  inches  was  produced  in  this  manner  and  stored  in  cold  storage  until 
ready  to  be  fabricated  into  specimens. 

The  procedure  for  layup  of  the  specimens  began  with  removing  the  prepreg 
from  cold  storage  and  cutting  it  to  proper  size  depending  on  the  type  of 
specimen.  In  all  cases,  a  dry  layer  of  181-style  fiber  glass  cloth  was 
used  as  a  bleeder  to  improve  resin- content  uniformity  and  to  insure  that 
uniform  venting  occurred  over  the  entire  laminate  surface  during  auto¬ 
clave  vacuum  and  pressure  application.  Teflon  FEP  film,  one-half  mil 
thick  and  perforated  with  3/ 16 -inch -diameter  holes  spaced  1/2  inch  apart 
in  a  rectangular  array,  was  used  as  the  parting  agent  between  the  bleeder 
and  the  prepreg. 

Regardless  of  the  type  of  specimen  (biaxial,  uniaxial,  or  tube),  a  B-stage 
time  of  approximately  7  hours,  from  time  of  removal  from  cold  storage  to 
time  of  placement  in  the  autoclave  for  curing,  was  used.  The  cure  used 
was  the  same  as  that  used  for  the  large- shell  facings  reported  in  Refer¬ 
ence  9;  namely,  160°  F  and  a  total  pressure  of  70  psi  CIO  psi  vacuum 
plus  60  psi  autoclave)  for  100  minutes,  with  a  zero-pressure  precure  of 
23  minutes.  Curing  was  conducted  in  the  large  autoclave  described  in 
References  9  and  11.  The  post  cure  was  350°  F  for  2  hours. 

To  fabricate  the  biaxial- load  specimens,  the  prepreg  was  arranged  in  the 
proper  lamination  arrangement  and  then  cut  to  conform  to  a  sheet-metal 
template.  The  template  has  a  small  bevel  cut  at  one  comer;  this  served 
to  index  the  prepreg  layers  during  final  layup  for  curing.  For  prepreg 
which  was  to  become  the  top  and  bottom  reinforcing  material  (i.e.,  all  of 
the  specimen  except  the  reduced- thickness  test  section),  a  circular  or 
elliptical  hole  was  cut  out  of  the  center,  the  shape  of  the  hole  depending 
upon  whether  or  not  the  specimen  was  to  be  for  a  1:1  or  a  1:2  biaxial-load 
ratio.  For  prepreg  which  was  to  become  the  test  section  and  the  central 
layers  of  the  rest  of  the  specimen,  no  hole  was  cut  out.  Special  slightly 
undersize  sheet-metal  inserts  which  fitted  into  the  cut-out  hole  were  used 
to  prevent  excess  resin  flow  into  the  reduced- thickness  test  section. 

Final  curing  in  the  autoclave  resulted  in  a  specimen  having  the  desired 
reduced- thicki jess  test  section  with  naturally  formed  fillets  of  resin. 
Next,  the  fillets  at  the  exposed  edges  between  the  loading  tabs  were 
machined;  and,  finally,  eight  loading  tabs  were  bonded  to  each  specimen. 
Bonding  was  accomplished  with  828-epoxy  and  V140-polyamide  curing  agent. 
The  bonding  cure  served  as  the  post  cure  for  the  torsion-tube  specimens. 
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The  uniaxial  tensile  specimens  were  cut  to  approximate  dimensions  from 
a  large,  completely  cured  sheet  of  the  proper  lamination  arrangement. 
Then  they  were  machined  to  final  dimensions. 

b.  Test  Procedures 


During  all  tests,  both  wet  and  dry  bulb  temperatures  were  measured. 

Also,  an  attempt  was  made  to  run  all  of  the  tests  at  approximately  the 
same  rate  of  loading. 

To  install  a  biaxial- load  specimen  in  the  loading  fixture,  a  portable 
hoist  was  used.  The  weight  of  the  specimen  plus  the  heavy  loading  tabs 
was  balanced  out  to  prevent  loading  the  top  part  of  the  specimen  in 
tension  and  the  bottom  part  in  compression.  Then  the  slack  was  pulled 
out  of  the  cable,  and  the  specimen  was  loaded  at  an  approximately  con¬ 
stant  rate.  Strain  readings  frcm  the  two  strain  gages  were  taken  on  a 
digital-diil  strain  indicator  at  predetermined  load  intervals. 

After  testing,  the  central  portion  of  the  specimen  was  cut  out  and  was 
saved  for  future  study.  Then  the  specimen  material  remaining  on  the 
loading  tabs  was  removed  by  sand-blasting  so  that  the  loading  tabs  would 
be  ready  to  be  used  again. 

Testing  of  the  tors ion- tube  and  uniaxial  tension  specimens  was  standard. 
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7.  EXPERIMENTAL  RESULTS 


This  section  presents  the  reduced  data  from  the  main  data- gathering  portion 
of  the  experimental  program.  The  experimental  stress-strain  or  load-strain 
data,  from  which  the  reduced  data  were  derived,  are  presented  in  Appendix  II. 

It  must  be  emphasized  that  the  stress  values  used  in  specifying  the  strengths 
are  nominal  composite  stresses,  i.e.,  total  load  divided  by  cross-sectional 
area.  These  are  the  values  that  an  aircraft  designer  needs  directly  in  design. 
However,  in  general,  they  are  not  the  same  as  the  stresses  physically  present 
in  each  individual  ply  of  the  laminate. 

a.  Uniaxial  Tensile  Tests 


The  individual  uniaxial  tensile  specimens  used  in  the  main  data- gathering 
series  were  designated  U1  through  U20.  Their  dimensions  are  shown  in 
Figure  10b.  The  data  reduced  for  individual  specimen  were  as  follows: 


Initial  Young's  modulus,  E 


Uniaxial  stress _ 

Strain  In  direction  of  stress 


(57) 


Initial  Poisson's  ratio,  v 


Transverse  strain _ 

Strain  in  direction  of  stress 


(58) 


Limit  strength  in  tension,  S,  =  Stress  at  0.01  percent  offset  strain 

L  (59) 

Ultimate  tensile  strength,  S,,  =  Stress  at  maximum  load  reached 

in  test  (60) 

The  stress  values  used  are  the  engineering  ones,  i.e.,  load  divided  by 
original  cross-sectional  area.  The  strain  values  used  are  those  measured 
by  variable-electric- resistance  strain  gages. 

The  initial  Young's  moduli,  E,  and  initial  Poisson's  ratios,  v,  were 
both  determined  in  the  initial  straight-line  elastic  portion  of  the 
stress-strain  curves. 

In  aircraft  metallic  alloys,  the  criterion  used  to  define  limit  strength 
is  the  0. 2-percent-offset  yield  strength.  However,  many  composite 
materials,  including  glass  fiber-reinforced  plastics  (FRP),  do  not  reach 
this  high  of  an  offset  strain.  Therefore,  for  such  materials,  a  new 
limit-strength  criterion  must  be  used.  The  limit- strength  criterion  used 
here  is  the  0.01-percent- off  set  yield  strength,  which  is  the  strain 
associated  with  an  offset  of  0.01  percent  strain. 

The  ultimate  tensile  strength  is  defined  in  the  traditional  manner.  For 
many  compos ite  materials,  including  FRP,  there  is  no  drop  in  load  with 
continuing  strains;  i.e.,  the  fracture  strength  coincides  with  the 
ultimate  tensile  strength. 
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The  reduced  data  for  the  main  series  of  uniaxial  tensile  tests  are  pre¬ 
sented  in  Table  TV.  (Additional  exploratory  uniaxial  data  are  presented 
in  Section  4c.)  It  is  noted  that  the  cross-p]y  specimens  were  not 
tested  at  a  90-degrees  orientation,  since  this  is  no  different  than 
0-degree  orientation  for  such  a  specimen  (symmetrically  laminated  with 
an  equal  number  of  plies  at  0  and  90  degrees) . 

Mean  values  of  the  moduli  (E  and  v) ,  used  in  determining  biaxial  stress 
values,  are  presented  in  Table  V. 

b.  Torsion-Tube  Shear  Tests 


The  individual  torsion-tube  shear  specimens  used  in  the  main  data- 
gathering  series  (see  Figure  10a)  were  designated  SI  through  S24.  The 
data  reduced  for  each  specimen  were  as  follows: 

Initial  shear  modulus,  G  =  shear  stress/shear  strain  (61) 

Limit  shear  strength,  S.,  =  shear  stress  at  0.01  percent 

offset  shear  strain  (62) 

Ultimate  shear  strength,  =  shear  stress  at  maximum 

torque  reached  in  test  (63) 

Since  the  lamination  arrangements  were  all  symmetric  about  the  middle 
of  the  thickness  and  since  the  specimens  were  thin-walled  and  circular 
in  cross  section,  the  following  expression  for  shear  stress  is  valid 
throughout  the  test: 

a6  =  TRq/J  =  2TRq/(R£  -  R$)  (64) 

where  T  =  torque,  J  =  polar  area  moment  of  inertia,  R  -  outside  radius, 
and  R^  =  inside  radius. 

For  a  state  of  pure  shear  strain,  application  of  the  Mohr  strain  circle 
shows  that  the  shear  strain  is  equal  to  twice  the  reading  of  an  electric- 
resistance  strain  gage,  which  measures  noxmal  strain  only,  at  an  orienta¬ 
tion  of  45  degrees  with  the  plane  on  which  the  shear  strain  takes  place. 
To  average  out  any  small  amount  of  bending  which  may  have  been  induced, 
two  gages  were  used,  one  at  +45  and  the  other  at  -45  degrees.  Thus,  the 
shear  strain  plotted  in  Appendix  I  is  merely  the  sun  of  the  strain  read¬ 
ings  from  the  two  gages. 

The  shear  moduli  and  shear  strengths  are  defined  in  a  fashion  analogous 
to  that  used  for  the  corresponding  uniaxial  quantities.  Figure  18  shows 
a  typical  shear  specimen  fracture. 

The  reduced  data  are  presented  in  Table  VI.  Mean  values  of  the  shear 
moduli  for  each  laminate  at  0-degrec  orientation  are  listed  in  Table  V. 
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TABLE  IV.  RESULTS  OF  UNIAXIAL  TENSILE  TESTS 


Specimen 

No. 

E,  psi 

V 

SL,  psi 

Su,  psi 

Four-Ply  Parallel  Laminated, 

0-deg  Orientation 

U1 

3.34  x  106 

0.148 

13,250 

37,400 

U2 

3.24  x  106 

0.144 

11,750 

41,150 

U3 

3.36  x  106 

0.132 

13,500 

40,500 

Four-Ply  Parallel  Laminated, 

90-deg  Orientation 

U4 

3.14  x  106 

0.193 

10,000 

41,400 

U5 

3.19  x  106 

0.181 

12,000 

42,600 

U6 

3.11  x  10u 

0.192 

12,750 

42,200 

Four-Ply  Cross  Laminated,  0-deg  Orientation 

U7 

4.55  x  106 

0.252 

11,750 

23,200 

U8 

3.10  x  106 

0.152 

11,250 

25,500 

U9 

3.16  x  106 

0.148 

14,750 

22,200 

U10 

3.33  x  106 

0.095 

10,500 

25,600 

Six-Ply  Quasi- Isotropic  Laminated,  0-deg  Orientation 

Ull 

3.59  x  10b 

0.297 

13,250 

33,700 

U12 

3.42  x  106 

0.324 

12,500 

35,000 

U13 

3.68  x  106 

0.320 

15,750 

31,100 

U14 

3.79  x  106 

0.323 

19,000 

34,000 

U15 

3.63  x  106 

0.323 

12,250 

39,600 

Six-Ply  Quasi -Isotropic  Laminated,  90-deg  Orientation 

U16 

3.63  x  106 

0.286 

10,500 

36,400 

U17 

3.85  x  106 

0.336 

13,250 

33,800 

U18 

3.62  x  106 

0.312 

13,500 

33,700 

U19 

3.57  x  106 

0.328 

12,000 

35,500 

’J20 

3.61  x  106 

0.291 

11,250 

34,500 
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TABLE  V.  SUMMARY  OF  MEAN  VALUES  OF  INITIAL  COMPOSITE  MODULI 


Lamination  Arrangement 

Parallel -Ply 

Cross -Ply 

Quasi -Isotropic 

Major  elastic  modulus, 

611 

3.304  x  106 

3.196  x  106 

3.62  x  106** 

Minor  elastic  modulus, 

£22 

3.15  x  106 

3.196  x  106 

3.66  x  106** 

Major  Poisson's  ratio, 

C12 

0.1410 

0.1316 

0.317*** 

Minor  Poisson's  ratio, 

V21 

0.188* 

0.1316 

0.311*** 

Shear  modulus,  G 

1.325 

1.31  x  106 

1.93  x  106 

*  This  value  appears  to  be  erroneously  high,  so  a  value  of  0.1342,  cal¬ 
culated  from  the  reciprocal  relationship  (0.1410  x  3.15/3.304)^35 
used  in  all  subsequent  data  reduction. 

**  Since  these  values  are  not  significantly  different,  within  experi¬ 
mental  error,  a  single  value  of  3.64  x  10^  psi  is  used  in  all  sub¬ 
sequent  data  reduction. 

***  Since  these  values  are  not  significantly  different,  within  experi¬ 
mental  error,  a  single  value  of  0.314  is  used  in  all  subsequent 
data  reduction. 


TABLE  VI.  RESULTS  OF  TORSION-TUBE  SHEAR  TESTS 


Specimen 

No. 

G,  psi 

^6’  Psi 

SU6>  psi 

Four-Ply  Parallel  Laminated,  0-deg  Orientation 

SI 

1.47  x  106 

7,000 

30,100 

S2 

1.25  x  106 

7,200 

31,256 

S3 

1.53  x  106 

4,000 

79,862 

S4 

1.05  x  106 

8,500 

28,040 

Four - Pi y 

Parallel  Laminated,  45 -deg  Orientation 

S5 

3.68  x  106 

6,800 

57,400 

S6 

3.13  x  106 

5,500 

54,000 

S7 

7.7  x  106* 

7,200 

80,800* 

S8 

3.37  x  106 

9,000 

66,500 

Four-Ply  Cross  Laminated, 

0-deg  Orientation 

S9 

1.3  x  106 

4,800 

28,125 

S10 

1.43  x  106 

7,000 

21,350 

Sil 

0.96  x  106 

10,200 

24,766 

S12 

1.75  x  106 

5,000 

32,260 

S13 

1.11  x  106 

7,000 

25,620 

S14 

2.3  x  106* 

10,500 

32,238 

Four-Ply  Cross  Laminated, 

45-deg  Orientation 

SIS 

2.5  x  106 

10,200 

50,000 

S16 

3.34  x  106 

10,500 

46,875 

SI  7 

2.22  x  106 

7,000 

46,875 

TABLE  VI  -  C 

ontinued 

Specimen 

No. 

G,  psi 

is-  psi 

is-  psi 

Six-Ply  Quasi 

-Isotropic  Laminated,  0-deg  Orientation 

S18 

2.0  x  106 

9,800 

55,510 

S19 

2.22  x  10° 

11,000 

53,375 

S20 

1.82  x  106 

10,500 

55,510 

Six -Ply  Quasi 

-Isotropic  Laminated,  45-deg  Orientation 

!  S21 

1.82  x  106 

7,500 

26,218 

S22 

1.62  x  106 

14,800 

49,105 

S23 

1.82  x  106 

10,000 

49,319 

S24 

2.22  x  106 

16,400 

50,386 

*  Values  too  high;  not  used 

in  final  data 

reduction. 
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The  value  for  the  minor  Poisson's  ratio  for  the  parallel  ply  presented 
in  Table  V  appears  to  be  high,  so  a  calculated  value  was  used  for  all 
data  reduction.*  Also,  the  values  for  the  major  and  minor  elastic 
moduli  were  not  different,  within  experimental  error,  so  a  mean  value 
was  used  for  each  of  them.  This  mean  value  was  taken  for  the  quasi¬ 
isotropic  and  was  an  average  of  the  values.  The  minor  and  major  moduli 
for  the  parallel  and  cross  ply  were  close,  within  experimental  error; 
therefore,  the  data  were  not  averaged  for  the  data  reduction. 


Biaxial -Loading  Tests 


The  individual  biaxial- loading  specimens  used  in  the  main  data- gathering 
series  (see  Figure  9)  were  designated  B1  through  B54,  and  the  correspond¬ 
ing  load-strain  curves  are  presented  in  Appendix  II. 


The  composite  biaxial-stress  values  were  obtained  by  multiplying  the 
measured  load  value,  P,  by  the  following  factors: 


oJ/P  «  (qpCeJ/P)  +  (Q'2)(eJ/P), 


(65) 


°2/P  =  (QJ2)  CeJ/p)  +  (Q$2)C4/P), 


(66) 


where  directions  1  and  2  are  the  major  and  minor  material-symmetry  direc¬ 
tions  (see  Section  2),  ej/P  and  e^/P  are  taken  from  the  biaxial-load 
versus  strain  data  presented  in  Appendix  II,  and  v,~  and  \>21  are  compo¬ 
site  Poisson's  ratios  taken  fran  Table  V.  For  0  =  0  degrees,  the  primes 
can  be  removed  and 

$il  =  $11  =  $k  =  $12  =  V21  ^11^>  $22  =  $22  =  ^22^  (-67-) 


where  A  =  1  -  V21 

For  0  =  45  degrees. 

$11  =  $22  =  (1/4^  ^11  +  2v21  hi  +  ^22  +  (-68) 

$12  =  *-1//4^  ^11  +  2v21E11  +  ^22  '  4^ 

Table  VII  presents  the  results  obtained  by  applying  Equations  (65)  and 
(66)  to  the  load-strain  data  of  Appendix  II. 

The  limit  and  ultimate  strength  criteria  are  the  same  as  those  used  for 
the  uniaxial  tensile  specimens.  There  appears  to  be  some  theoretical 


*The  calculated  value  was  obtained  from  the  reciprocal  relationship: 

0.1410  X  3.15/3.304. 
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Table  vii  .  results  of  biaxial 

-LOADING  TESTS 

Specimen 

No. 

Sj^ ,  psi 

Syj ,  psi 

SL2*  psi 

SU2’  psi 

Four-Ply  Parallel  Laminated,  0-deg  Orientation,  1:2  Nominal  Loading 

B12 

18,900 

33,000 

9,700 

14,900 

B13 

19,650 

26,500 

9,150 

12,420 

B14 

15,800 

29,000 

6,370 

12,150 

B15 

17,900 

35,100 

7,560 

15,550 

B16 

14,980 

28,500 

6,840 

13,750 

Four-Ply  Parallel  Laminated,  0-deg  Orientation,  1:1  Nominal  Loading 

B1 

13,200 

29,500 

13,280 

28,400 

B2* 

18,000 

28,000 

17,200 

30,800 

B3* 

16,820 

30,850 

14,650 

29,650 

B4 

16,700 

28,450 

15,500 

29,250 

B5 

20,500 

28,200 

15,850 

29,000 

B6 

15,750 

36,000 

12,850 

37,000 

B17 

12,950 

33,700 

15,200 

35,900 

Four-Ply  Parallel  Laminated,  0-deg  Orientation,  2:1  Nominal  Loading 

B18 

6,800 

11,700 

12,300 

27,000 

B19 

7,650 

13,000 

25,700 

36,300 

B20 

10,350 

13,800 

18,000 

28,400 

B21 

8,320 

14,650 

19,300 

35,100 

B22 

11,300 

14,250 

20,900 

29,750 

Four-Ply  Parallel  Laminated,  45-deg  Orientation,  1:1  Nominal  Loading 

B7 

25,200 

38,800 

23,200 

38,800 

B8 

15,100 

24,200 

13,600 

28,400 

B9* 

10,950 

26,500 

14,350 

28,000 

BIO* 

14,100 

23,800 

11,040 

25,100 

Bll 

10,250 

20,850 

16,100 

29,300 
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TABLE  VII  -  Continued 

Specimen 

No. 

SL1>  Psi 

Sy^ ,  psi 

SL2’  Psi 

S(J2 ’  Psi 

Four-Ply  Cross 

Laminated , 

0-deg  Orientation 

,  1:1  Nominal  Loading 

B23* 

14,200 

23,000 

14,350 

26,960 

B24 

15,200 

31,700 

16,700 

31,700 

B25 

19,100 

26,500 

19,100 

26,500 

B26 

19,100 

27,200 

20,950 

27,200 

B27 

18,400 

26,400 

15,400 

26,400 

B28 

17,650 

28,400 

17,650 

28,400 

Four-Ply  Cross 

Laminated, 

0-deg  Orientation 

,  1:2  Nominal  Loading 

B29* 

6,800 

13,250 

16,700 

30,400 

B30 

9,520 

16,350 

21,000 

38,200 

B31 

11,400 

18,900 

13,300 

23,000 

B32* 

6,750 

10,200 

14,100 

19,400 

B33* 

4,000 

4,000 

9,000 

10,200 

B34 

9,370 

13,350 

16,000 

25,300 

B35* 

6,790 

11,430 

15,250 

26,700 

Six-Ply  Quasi- 

Isotropic  Laminated,  0-deg  Orientation, 
Loading 

1:1  Nominal 

B36 

12,600 

26,850 

13,580 

28,000 

B37* 

13,500 

39,100 

9,450 

32,500 

B38 

10,650 

34,100 

11,700 

34,100 

B39 

13,000 

38,600 

16,200 

41,600 

B40 

5,930 

24,800 

6,000 

26,900 

Six-Ply  Quasi- 

Isotropic  Laminated,  45-deg  Orientation, 
Loading 

1:1  Nominal 

B41* 

14,200 

27,700 

14,950 

27,700 

B42 

7,300 

14,000 

11,850 

19,500 

B43 

15,600 

28,300 

17,350 

24,700 
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TABLE  VTI  -  Continued 


Specimen 

No. 

SLi’  Psi 

S^,  psi 

SL2’  Psi 

y,2,  psi 

B44 

12,350 

20,800 

14,300 

22,200 

B45 

18,950 

13,200 

20,900 

Six- Ply  Quasi- 

Isotropic  Laminated,  0-deg  Orientation, 
Loading 

2:1 

Nominal 

B50* 

20,700 

33,900 

10,600 

13,350  1 

B51 

25,300 

36,800 

13,050 

13,900  : 

B52* 

22,000 

41,200 

12,200 

16,600 

BS3* 

21,100 

34,600 

11,350 

12,400 

B54* 

22,700 

41,000 

14,600 

15,300 

Six-Ply  Quasi- 

Isotropic  Laminated,  0-deg  Orientation, 
Loading 

1:2 

Nominal 

B46* 

7,150 

12,600 

21,950 

36,000 

B47 

7,360 

9,700 

18,600 

23,600 

B48* 

6,600 

11,350 

14,500 

33,000 

B49* 

9,600 

13,420 

20,600 

36,000 

*  In  each  specimen  marked  with  an  asterisk 
pass  through  the  reduced -thickness  test 

,  the  fracture  did  not 
section. 
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basis  for  using  an  "equivalent"  offset  strain  for  biaxial  loading.  How¬ 
ever,  practical  design  usage  appears  to  favor  using  the  same  offset  for 
all  types  of  loading,  as  discussed  by  Bert  et  al  (Reference  36). 

It  should  be  mentioned  that  the  method  used  here  to  determine  the  stress 
from  measured  strains  is  quite  accurate  in  the  linear  elastic  range. 

When  the  load-strain  curve  becomes  nonlinear,  the  accuracy  of  the  stress 
values  is  not  known.  There  are  two  factors,  both  of  which  would  tend 
to  reduce  the  actual  stress  to  a  lower  value  than  predicted  by  the  calcu 
lation.  If  the  material  in  the  reduced  test  section  is  yielding  (i.e., 
undergoing  higher  strains  than  it  would  if  it  were  elastic) ,  its  effec¬ 
tive  moduli  would  be  lower  than  the  corresponding  elastic  moduli  and 
thus  would  lower  the  actual  stress.  However,  if  material  not  in  the 
reduced- thi ckness  test  section  yields  first,  additional  load  may  be 
dumped  into  the  test  section,  causing  the  strain  gages  to  read  errone¬ 
ously  high  for  a  given  value  of  load  applied  to  the  edge  of  the  specimen 
Thus,  the  calculated  ultimate  strength  values  presented  represent  upper 
bounds. 

Figure  19  shows  photographs  of  typical  fracture  patterns  for  both  speci¬ 
men  types  (1:1  and  1:2  or  2:1).  The  fact  that  a  number  of  the  fractures 
of  the  1:1  type  specimens  were  multiple  is  additional  testimony  to  the 
very  good  attainment  of  a  uniform  balanced-biaxial  stress  field  in  the 
test  section.  Of  the  54  biaxial- loading  specimens  tested  in  the  data- 
gathering  series,  only  18  (marked  with  an  asterisk  in  Table  VII)  did  not 
fracture  through  the  test  section.  However,  a  study  of  the  apparent 
ultimate  strength  data  in  Table  VII  shows  that  these  specimens  failed 
at  values  very  close  to  those  of  the  other  specimens,  i.e.,  neither 
abnormally  high  nor  abnormally  low.  This  indicates  that,  in  most  cases, 
test-section  failure  was  probably  imminent. 
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Figure  19.  Photographs  of  Typical  Fracture  Patterns  for  Biaxial -Loading 
Specimens . 
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EVALUATION  OF  EXPERIMENTAL  RESULTS 


In  this  section,  the  experimental  results  presented  in  Section  7  are  evalu¬ 
ated  in  the  light  of  the  theoretical  stiffness  and  strength  considerations 
presented  in  Sections  2  and  3. 

a •  Parallel-Ply  Laminates 

Since  there  is  presently  no  composite  micromechanics  analysis  known  to 
the  authors  to  be  applicable  to  fabric-reinforced  composites,  no  quanti¬ 
tative  prediction  of  the  individual  stiffnesses  of  parallel -ply  laminates 
can  be  made  on  the  basis  of  the  properties  of  the  constituent  materials 
(glass  and  epoxy,  in  this  case.)  However,  since  the  reciprocal  relation¬ 
ship  is  based  on  very  fundamental  thermodynamics  principles,  as  discussed 
in  Section  2a,  it  should  hold.  The  two  parallel -ply  moduli  values  pre¬ 
sented  in  Table  V  appear  to  be  reasonable  for  the  resin  content  (34-36 
percent  by  weight;,  and  so  does  the  value  of  v^.  However,  the  \>2-  value 
appears  to  be  erroneous,  and  since  the  reciprocal  relation  was  violated, 
a  value  of  0.1342  given  bv  the  reciprocal  relation  appears  to  be  more 
valid.  The  v^2  and  the  corrected  as  well  as  the  E,,  and  E™  values, 
are  quite  close  to  those  reported  by  Youngs  (Reference  80)  and  Military 
Handbook  17  (Reference  81) . 

The  value  of  G  presented  in  Table  V,  1,386,000  psi,  is  approximately  71 
percent  higher  than  the  810,000-psi  value  given  in  Military  Handbook  17. 
Part  of  the  difference  may  be  due  to  the  type  of  specimen,  since  the 
Military  Handbook  17  values  were  obtained  by  in-plane  panel-shear  tests. 
Sidorin  (Reference  82)  obtained  11  percent  higher  shear-modulus  values 
by  torsion-tube  tests  than  by  panel-shear  tests  for  the  same  fabric 
orientation  under  discussion  here. 

The  value  of  5'  corresponding  to  an  orientation,  0,  of  45  degrees  can  be 
predicted  from  the  values  of  5,  E,-.,  E??,  and  v,_,  by  adapting  the  last 
of  Equations  (17)  as  follows: 

1/G'  =  (1/5)  +  [(1  +  2v12)  ExJ  +  (1/E22)  -  (1/5)]  sin2  2e  (69) 

Equation  (69)  gives  a  value  of  1,37Q,00Q  psi,  which  is  considerably 
lower  than  the  3,39Q,QQQ-psi  value  of  Er  determined  experimentally  for 
0  =  45  degrees.* 

Figure  20  presents  the  limit  and  apparent  ultimate  strength  results  from 
Table  VII  in  the  form  of  strength  envelopes.  It  is  noted  that  there 


"It  appears  that  the  incompatibility  of  the  torsion-tube  properties  and  the 
properties  of  the  flat  specimens  (uniaxial  and  biaxial)  was  due  to  a  differ¬ 
ence  in  resin  content.  The  tors ion- tube  specimens  did  not  bleed  adequately, 
resulting  in  a  resin  content  of  approximately  40  percent,  compared  with  34 
percent  for  the  flat  specimens. 
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appears  to  be  considerable  symmetry  about  the  1:1  biaxial  data;  i.e., 
the  2:1  data  and  1:0  data  arc  symmetric  with  the  1:2  and  0:1  data, 
respectively.  No  attempt  was  made  to  fit  a  strength  theory  to  the 
data  presented. 

b.  Cross -Ply  Laminates 

In  Section  2d,  it  was  shown  that  laminate  stiffness  theory  predicts  that 
the  stiffness  properties  of  a  symmetrically  laminated  cross-ply  laminate 
are  simply  the  mean  values  of  the  properties  of  the  individual  plies. 

Thus,  the  predicted  values  would  be  as  follows: 

i;n  =  U7?  =  3.23  x  106  psi,  v.7  =  \>22  =  0.1376 

Both  of  these  values  are  slightly  higher  than  the  corresponding  values 
listed  in  Table  V  for  the  cross-ply  laminate. 

Figure  21  shows  the  limit  and  apparent  ultimate  strength  envelopes  for 
the  cross-ply  laminate,  drawn  from  data  in  Table  VII. 

c.  Quasi -Isotropic  Laminates 

'Hie  predicted  quasi -isotropic  stiffness  values  obtained  by  using  parallel- 
ply  stiffness  properties  in  liquations  (37),  Section  2,  are  as  follows: 

Eh  =  =  3.04  x  106  psi, 

v12  =  V21  =  °’160’ 

C  =  1.31  x  106  psi. 

These  values  compare  very  poorly  with  the  mean  quasi -isotropic  values 
listed  in  the  last  column  of  Table  V. 

Figure  22  shows  the  limit  and  apparent  ultimate  strength  envelopes  for 
the  quasi -isotropic  laminate,  using  values  from  Table  VII.  The  obviously 
higher  yield  strength  at  a  biaxial-load  ratio  of  2:1  as  compared  to  1:1 
is  not  explained. 
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Figure  22.  Biaxial -Tension  Strength  Envelope  for  6-Ply  Quasi- Isotropic 
Laminate. 


9.  CONCLUSIONS  AND  RECOMMENDATIONS 


a.  Conclusions 

(1)  Use  of  the  biaxial- loading  specimens  and  associated  biaxial- loading 
fixture  developed  in  this  program  is  believed  to  be  the  most  accur¬ 
ate,  reliable,  and  simple- to- fabricate  method  for  obtaining  strength 
data  on  flat-sheet  composite  materials  subjected  to  biaxial- tens  ion 
loading . 

(2)  The  Volan-A-finish,  E-glass,  181-style  fabric/828-Z  eopxy  laminates 
tested  in  this  program  generally  exhibited  considerably  nigher  limit 
strengths  (and  slightly  lower  apparent  ultimate  strengths)  under 
1:1,  1:2,  and  2:1  biaxial-load  ratios  than  in  uniaxial  tension. 

(3)  Of  the  three  lamination  arrangements  evaluated  (parallel-ply, 
cross-ply,  and  quasi- isotropic) ,  the  parallel-ply  laminate  was 
slightly  superior  under  biaxial  loading,  although  the  differences 
were  small.  This  suggests  that  for  this  particular  weave,  combina¬ 
tion  of  materials,  and  type  of  loading  (uniaxial  or  biaxial  tension), 
parallel-ply  lamination  is  advantageous  structurally  as  well  as 
being  easier  to  fabricate.  (This  same  conclusion  would  not  be 
expected  to  hold  true  for  a  more  highly  ’mdirectional  composite 

such  as  a  tape  or  a  filament- wound  composite.) 

(4)  The  thin-walled  torsion-tube  specimen  used  in  this  program  is 
believed  to  be  an  accurate  and  reliable  method  for  obtaining  shear 
modulus  and  shear  strength  data  on  composite  materials.  However, 
it  does  suffer  from  the  requirement  that  the  composite  must  be 
fabricated  into  the  form  of  a  tube. 

b .  Re  commendat  i  ons 


(1)  Additional  composite  materials  of  importance  for  Army  aircraft 
primary  structures  should  be  evaluated  in  biaxial-tension  loading 
using  the  types  of  specimens  and  loading  device  developed  in  the 
program.  In  view  of  Conclusion  (3),  it  would  be  especially  impor¬ 
tant  to  evaluate  a  more  highly  unidirectional  material  such  as  a 
tape.  The  fiber  material  may  be  commercial  S-glass,  boron,  silicon 
carbide,  graphite,  etc. 

(2)  In  conjunction  with  Recommendation  (1) ,  it  is  recommended  that  the 
optimal  lamination  arrangement  be  determined  for  each  combination 
of  composite  material  and  biaxial- loading  ratio. 

(3)  In  view  of  the  strong  dependence  of  uniaxial  tensile  strength  values 
on  specimen  width,  length,  number  of  plies,  and  presence  of  stress- 
concentration-producing  discontinuities,  such  as  changes  in  thick¬ 
ness,  it  is  recommended  that  more  research  be  devoted  to  this 
problem.  A  more  extensive  experimental  program,  coupled  with  an 
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appropriate  analytical  investigation,  appears  to  be  advisable  in 
order  to  arrive  at  sufficiently  definitive  information  for  more 
reliable  design  of  uniaxially  loaded  structures  of  composite 
materials. 
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APPENDIX  I 


DESIGN  OF  TEST-SECTION  CONiOUR  TO  ACHIEVE 
A  UNIFORM  STRESS  ‘'STa’IT W ffi'STREt)"BOTTAL'-,LOAl)  RATIO 


As  pointed  out  in  Section  4  of  the  body  of  the  report,  for  material-property 
evaluation,  it  is  desirable  to  have  a  uniform  stress  state  in  the  test  sec¬ 
tion.  In  this  appendix,  the  shape  of  reduced- thickness  test  section  neces¬ 
sary  to  achieve  a  uniform-biaxial-stress  state  is  derived  by  using  two- 
dimensional  theory  of  elasticity. 

1 .  Hypotheses 

The  following  hypotheses  are  made  in  the  ensuing  synthesis: 

CHI)  Generalized  plane  stress  conditions  exist  throughout  the  reduced¬ 
thickness  test  section  and  the  reinforced  section;  i.e.,  all  stress 
components  in  the  thickness  direction  are  neglected. 

(H2)  The  only  loadings  considered  are  edge  loadings;  i.e.,  all  body 
forces  such  as  those  due  to  gravity  and  inertia  are  neglected. 


2.  Synthesis 

To  assure  that  the  equations  of  equilibrium  of  a  body  under  generalized- 
plane-stress  conditions  are  satisfied  identically,  it  is  customary  to 
introduce  the  Airy  stress  function,  <t> ,  defined  as  follows: 


(70) 


where  N  ,  N  ,  N  are  the  internal  forces  per  unit  length  of  run  (normal 
to  the  $  an&  y  Sees,  and  shear,  respectively),  and  a  subscript  comma 
denotes  partial  derivatives  with  respect  to  the  variables  appearing  in 
the  subscript  after  the  corrcna. 


A  uniform  stress  state  requires  that  N  ,  N  ,  and  N  all  be  constants. 
Then  integration  of  Equations  (70)  yields  ¥he  following  expression  for 
the  Airy  stress  function  associated  with  such  a  stress  state: 

t  =  (1/2)  (N  y2  +  N  x2)  -  N  xy  +  Ax  +  By  +  C  (71) 

a  y  a/ 

where  x  and  y  are  Cartesian  position  coordinates,  and  A,  B,  and  C  are 
arbitrary  constants  of  integration. 


The  required  contour  shape  (position  coordinates  xc,  yc)  of  reduced¬ 
thickness  test  section  is  the  one  in  which  the  stress  function  is 
invariant  with  respect  to  position  xc,  y  .  This  is  satisfied  by  the 
following  expression: 

-N  xy+Ax+By+K  =  0 
xy  e'e  c  7c 


(1/2) (Nx  y2  +  Ny  x2) 
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(72) 


where  K  =  <f>(*c,  yc)  -  C. 

The  constants  A  and  B  in  Equation  (72)  can  be  omitted,  since  they  merely 
determine  the  origin  of  the  coordinate  system. 

Furthermore,  for  any  generalized-plane- stress  state  (Nx,  I\L,  N^),  by 
means  of  the  Mohr's  stress  circle,  an  orientation  can  be  found  such  that 
the  shear  force  N  vanishes.  The  loadings  associated  with  this  orien¬ 
tation  are  known  ns  the  principal  loadings  N_,  N„.  If  the  axes  of  the 
Cartesian  coordinate  system  used  to  define  the  shape  of  the  reduced¬ 
thickness  test  section  are  chosen  to  coincide  with  the  principal  loading 
directions,  can  be  omitted  and  Nx  and  N  are  replaced  by  NL  and  Nq. 
Thus,  Equation  (72)  can  he  rewritten  as  follows: 

(1/2)  (Np  V2C  +  Nq  x2)  +  K  =  0  (73) 

Equation  (73)  can  be  rewritten  in  the  standard  fom  for  an  ellipse  as 
follows: 


(xc/a)2  +  (yc/b)2  =  1  (74) 

where  a  and  b  are  the  major  and  minor  semi -axes,  respectively. 

Comparing  Equations  (73)  and  (74) ,  it  is  found  that  the  following 
relationship  must  hold: 

b/a  =  (Nq/Np)1/2  (75) 

Fur  the  mo  re,  K  is  a  scale  factor  given  by 

K  =  -  N  N  /2  (76) 

P  q 

Applying  Equation  (75),  it  is  seen  that  for  a  balanced  biaxial  loading 
(biaxial- load  ratio  Nq/Np  =  1),  b/a  =  1;  i.e.,  the  required  contour  of 
the  reduced  thickness  section  is  a  circle.  For  a  biaxial- load  ratio  of 
1/2,  b/a  =  (1/2) 1/2  =  0.707. 

It  is  important  to  note  that  only  equilibrium  considerations  have  been 
used;  thus,  the  results  are  completely  independent  of  the  nature  of 
the  elastic  behavior  (anisotropic,  orthotropic,  or  isotropic). 
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Figure  25.  Stress- Strain  Plot  of  Uniaxial,  Cross  4-Ply,  0°  Orientation. 
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Load-Strain  Data  from  Biaxial- Loading  Tests 


t - 1 - 1 - 1 - ' - r 


iA  O  i/>  O  »A  O 

IN  CM 


Ol  x  ai  ~  oven 


97 


98 


Figure  36.  Load-! 
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Figure  37.  Load- 
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Figure  39.  Load-Strain  Plot  of  Biaxial,  1:2,  Parallel  4 -Ply,  0°  Orientation. 
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